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Abstract. Experimental form factors of the hydrogen and helium isotopes, extracted 
from an up-to-date global analysis of cross sections and polarization observables 
measured in elastic electron scattering from these systems, are compared to predictions 
obtained in three different theoretical approaches: the first is based on realistic 
interactions and currents, including relativistic corrections (labeled as the conventional 
approach); the second relies on a chiral effective field theory description of the strong 
and electromagnetic interactions in nuclei (labeled yEFT); the third utilizes a fully 
relativistic treatment of nuclear dynamics as implemented in the covariant spectator 
theory (labeled CST). For momentum transfers below Q < 5 fm“^ there is satisfactory 
agreement between experimental data and theoretical results in all three approaches. 
However, at Q > 5 fm“^, particularly in the case of the deuteron, a relativistic 
treatment of the dynamics, as is done in the CST, is necessary. The experimental 
data on the deuteron A structure function extend to Q ~ 12 fm“^, and the close 
agreement between these data and the CST results suggests that, even in this extreme 
kinematical regime, there is no evidence for new effects coming from quark and gluon 
degrees of freedom at short distances. 
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1. Introduction 

Few-nucleon systems, and more generally light s- and p-shell nuclei with mass number 
up to d = 12, offer a unique opportunity to test our understanding of nuclear dynamics. 
Over the past thirty years or so, several techniques have been developed to solve 
exactly the quantum mechanical few-body problem in both non-relativistic OEl and 
relativistic PHE] regimes. More recently, in the last decade, quantum Monte Carlo 
methods coupled with improvements in algorithms and advances in computational 
capabilities have made it possible to carry out also exact, albeit non-relativistic, 
calculations of nuclei up to PIT]. These technical breakthroughs have permitted 
hrst-principle studies of the strong interaction in nuclei, as it manifests itself in terms of 
two- and many-body forces among the nuclear constituents, the protons and neutrons, 
and of the interactions of these constituents with external electroweak probes in a wide 
range of energy and momentum transfers. 

In the present review the focus is on the electromagnetic ground-state structure of 
the hydrogen and helium isotopes. Since the early hfties, the associated form factors 
have been the subject of intense experimental and theoretical scrutiny. The large body 
of elastic electron scattering cross section (and polarization) data from these systems— 
a review of these data, and ensuing analysis, is provided in section |2]— has now led to 
accurate experimental determinations of the charge and magnetic form factors of ^H, 
^H, and ^He, the quadrupole form factor of ^H, and the charge form factor of '^He, up 
to momentum transfers Q beyond 9 fm“^, and in some instances, as for the A structure 
function of the deuteron, extending to Q ~ 12 fm“^. 

A vast amount of theoretical work in a variety of different frameworks—purely 
non-relativistic, or including relativistic corrections, or fully covariant ones—exists for 
these systems, and no attempt will be made here to summarize it. Rather, we have 
chosen to focus on three representative approaches: the hrst one, which we label as 
“conventional”, is based on realistic nuclear interactions and currents, including leading 
relativistic corrections; the second one relies on chiral effective held theory (yEFT) for a 
description of the nuclear strong and electromagnetic interactions; the third one utilizes 
the fully relativistic dynamical framework of the covariant spectator theory (CST). They 
are reviewed in considerable detail in section |3l where an appraisal of their differences 
and similarities is also provided. All three approaches have recently been used to 
calculate the few-nucleon form factors: conventional and yEFT predictions are available 
for A = 2-4 and CST ones for A = 2 and 3. These predictions are compared to the 
experimental form factors extracted from the world-data analysis of cross sections and 
polarizations in section 01 The yEFT calculation of the ^He form factor are presented 
here for the hrst time. Some concluding remarks are given in section [5l For completeness, 
in the remainder of this section we recall the basic formalism used to describe elastic 
electron-nucleus scattering and the dehnitions of the various form factors. 
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In the one-photon-exchange approximation, the nnpolarized cross section for elastic 
electron-deuteron scattering can be written as 

^ = f~l [A{Q) + B{Q) tan^ 0/2] = ctm fr~ell{Q, 0) (1-1) 


where ctm is the Mott cross section, 


(Tm — 


a cos 0/2 
_2e sin2 0/2 


( 1 . 2 ) 


/j.g/ is the recoil factor, 

/rec = 1 +sin^ 0/2 , (1.3) 

a is the fine structure constant, e and 0 are respectively the initial electron energy and 
hnal electron scattering angle, and Md is the rest mass of the deuteron. The structure 
functions A{Q) and B{Q) are functions only of the four-momentum transfer Q defined 
as Q = \/^ and = 4e£'sin^0/2, where e' is the final electron energy, and in elastic 
scattering the electron energy transfer e — e' is related to via e — e' = Q‘^/{2 Md). For 
a spin 1 nucleus like the deuteron, these structure functions can be expressed in terms 
of the three form factors Gc{Q), Gm{Q), and Gq{Q) (respectively charge, magnetic, 
and quadrupole form factor) as 


A{Q) — Gq{Q) + 2 Gm(Q) + 9 > (1-4) 

and 

BiQ) = +Vd)Gl^{Q) , (1.5) 

where rjd = Q^/(4Mj). Note that while Gm{Q) is given uniquely by B{Q), Gc{Q) 
and Gq{Q) both appear in A{Q) and, therefore, cannot be separated in an nnpolarized 
scattering experiment. Separation of the charge and quadrupole form factors requires 
elastic scattering involving polarization of either the initial or final deuteron states. 
For a tensor polarized initial deuteron, the tensor analyzing power T 2 o{Q) can be 
experimentally determined via 


I{Q,0)T2o{Q) 


^r^VdGc{Q)GQ{Q) + -r,^,Gl{Q) 

+ 2 + 2 (1 -|- r]d) tan^ 0/2] G\^{Q) 


( 1 . 6 ) 


where I{Q,9) is dehned in equation fll.ll) . This observable is especially sensitive to the 
ratio Gq{Q)/G ciQ)- So far, the separation of the three independent form factors has 
been carried out experimentally up to a momentum transfer Q — 8 fm“^. However, 
data on A{Q) extend up to Q — 12 fm“^. 

The form factors defined above can be related to matrix elements of the 
electromagnetic current between initial and final deuteron states. Introducing 
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Jxy = where e\^ is the virtual photon polarization, we dehne these matrix 

elements as 


(P', A' I .h_, IP, A) = g^^GUP', P)eX = G^l^iP', P) 


(1.7) 


and I P, A) and | P\ \') are the initial and hnal deuteron states with four momenta P^ 
and P^' and helicities A and A', respectively. Since there are only three scalar form 
factors, only three of the helicity matrix elements are unique. These are chosen to be 

9-1 = -4y- {P\ -1 l^ol P. 1) = ii(P', P) , (1.8) 

.0 (1.9) 

««" i 

On the other hand, by invoking Lorentz invariance, parity conservation, and current 
conservation, the matrix elements Gy^{P',P) can be shown to have the general form 

[SI El CD]: 


g;,,(p',p) 


- Gi(g) Vy<x-G,{q S'f^l"^ 

-G2(Q)(SG-,-^r &.«) , 


(p/ ^ py 


( 1 . 11 ) 


where the four-momentum transfer is g = P' — P, and = ^x{P) [■Ca' = ^x'{P')] 
the initial (hnal) deuteron helicity four-vectors. It then follows that 


9-1 — x/^ + Vd Gi{Q) , (1.12) 

9o = x/^ + Vd +‘^VdfGii^Q) — 2r](iG2{Q) + 2r](i{l + r]ci) Gs{Q)] , (1-13) 

9+1 = x/vd{^ +Vd) G2{Q) . (1.14) 


The charge, magnetic, and quadrupole form factors are related to the invariant functions 
Gi(Q)’s by 

GciQ) =G^{Q) + ^VdGQ{Q) , (1.15) 

Gm{Q)=G2{Q) , (1.16) 

Gq{Q) = Gi(Q) - G2{Q) + (1 + Vd) GsiQ) , (1.17) 

or directly to the matrix elements gx's by 

GciQ) = ir-7f^= i9o +‘^ 9-i) ) ( 1 . 18 ) 

3 V1 + % 

Gm(Q) = ! ^=9+1 ) (1.19) 

Vhd(l + T]d) 

GciQ) =- - (90 - 9-1) , ( 1 . 20 ) 

2 rjdx/1 + Vd 

and are normalized to 

Gc(0) = 1 , Gm(0) = ^t^d, Gq(0) = Ml Qd , (1.21) 
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where m is the nucleon mass, and and Qd are respectively the magnetic moment (in 
units of nuclear magnetons) and the quadrupole moment. In some of the calculations 
carried out in the conventional approach discussed in section 01 boost corrections in the 
initial and hnal states are ignored, and hence the factor 1 /on the right-hand-side 
of equations (I1.18p - (ll.20p is not included. 

Electron elastic scattering cross sections corresponding to a spin 0 target, like ^He, 
or spin 1/2 target, like ^He/^H, are well known [11], and will not be given here. The 
former can be expressed in terms of a single charge form factor, while the latter involves 
a charge and a magnetic form factor. Ignoring again boost corrections, the “^He charge 
form factor is 

fc(Q) = l (0|p(q)|0) , (1,22) 

while the trinucleon charge and magnetic form factors follow from 

FciQ) =l(l/2, + |p(q)|l/2,+) = F,(Q)-^f,(Q), (1,23) 

2 Til 1 

Fm{Q) =- (1/2, + b.(q)| 1/2, -) = -[Fi((3) + F,_{Q)] , (1,24) 

P fl 

where |0) represents the '‘He ground state, |l/2, ±) represent the trinucleon ground states 
with spin projections ±1/2 along the direction of the momentum transfer q, and Fi and 
±2 are the three-body Dirac and Pauli form factors that appear in the relativistic current 
for a target of mass Mt- The time and space parts of the four current introduced 
earlier are denoted respectively as p and j. In equations fll.221) and fll.23p - fll.24p Z is 
the proton number, p the trinucleon magnetic moment in units of nuclear magnetons, 
and m is the nucleon mass. Hence, the Fc{Q) and Fm{Q) form factors are normalized 
to 

FciO) = Fm(0) = 1 . (1.25) 

Below we will also consider the isoscalar and isovector combinations of the trinucleon 
charge and magnetic form factors, defined as (suppressing the Q-dependence for 
simplicity) 

F^/^ = i[2Fa('He)±Fc(=‘H)] , (1.26) 

^ [/i('He)FM('He) ± pC11)FmCB)] . (1.27) 

If the ^H and ^He ground states were pure isospin T = 1/2 states, then Fq,F^ and 
Fq,F^ would only be affected by, respectively, the isoscalar (S') and isovector {V) 
components of the current. However, small isospin admixtures with T > 1/2, induced by 
the electromagnetic interactions as well as charge-symmetry breaking terms in the strong 
interactions, are included in trinucleon wave functions calculated in the conventional and 
yEFT approaches discussed below. As a consequence, isoscalar (isovector) currents give 
non-vanishing, albeit small, contributions to the isovector (isoscalar) form factors. 
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2. Overview of the world experimental data 
2.1. Determination of form factors 

For the light nuclei of interest to this review, an extensive set of data covering a large 
range of electron energies e and scattering angles 6 is available. Such a set for various 
reasons is not very useful for a direct comparison to theory. Several considerations need 
to be re-emphasized: (i) for most nuclei the cross sections depend on two form factors, 
related to the Coulomb monopole Cq {Gc for the deuteron or Fc for the three- and 
four-body nuclei) and magnetic dipole Mi {Gm or Fm)- For the case of the deuteron 
a third form factor, due to the Coulomb quadrupole operator C 2 (Gq), is contributing 
as well. As compared to the cross sections, the individual form factors are much more 
sensitive to the ingredients of the theoretical calculations. It is therefore highly desirable 
to extract from the experimental cross sections the individual form factors, (ii) For the 
A = 3 nuclei (J^ = 1/2), the charge and magnetic form factors have traditionally 
been determined by the authors who performed the individual experiments. For the 
deuteron case (Jj = 1) the data have mainly been discussed in terms of the A{Q) and 
B{Q) structure functions and the observable T 2 o(Q), dehned in equations (11.dh . (II.hh . 
and fll.bp . 

To get the most precise form factors from the world data, it is necessary to re-analyze 
all cross sections and polarization observables at a given momentum transfer, as only 
the combination provides an optimal separation of the form factors. Furthermore, it is 
only with such a separation using the world data that a reliable error bar of the form 
factors can be derived. 

Re-analysis of the world cross sections is required for another reason: the step from 
cross section to form factor should involve removal of the Coulomb distortion of the 
electron waves, so that the form factors can be compared to the one-photon exchange 
results from theory. This step, mainly important at low momentum transfer, has been 
omitted in most of the past analyses determining form factors. In the re-analysis of the 
world data the Coulomb corrections, which depend on both e and 6, can be accounted 
for. They are done using the approach described in Ref. [12]. 

A further complication arises from the fact that the cross sections and analyzing 
powers have been measured at a variety of energies and angles. In general even 
individual experiments aiming at a longitudinal/transverse (L/T) separation have not 
achieved exactly the same momentum transfers at forward and backward angles, so 
non-transparent interpolations/extrapolations are necessary; this difficulty is even more 
serious when combining data from different experiments. 

In order to get the most precise information, we proceed as follows: starting from 
the experimental world cross sections data, we correct them for Coulomb effects (and, 
if desired, the complete set of two-photon exchange corrections) and fit the resulting 
one-photon exchange cross sections using a very flexible parameterization for the two 
(three) form factors of interest. For any interval of momentum transfer, this provides 
the most complete information on both the separated form factors and their uncertainty. 
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In this fitting procedure, the statistical uncertainties of the data can be included 
in the standard way using standard error propagation and the error matrix. For the 
determination of the systematic errors, we employ a rather conservative approach: the 
data from individual experiments are changed by the quoted systematic error. Then 
the resulting changes of the £t form factors are quadratically added and quadratically 
combined with the statistical uncertainties. 

While the results from such a fit do provide the best experimental information, 
one peculiarity needs to be understood: the values of the form factors at closely spaced 
values of momentum transfer are not independent. When htting the cross sections with 
parameterized form factors, the correlation extends over an interval AQ ~ 1/Rmax, 
where Rmax is the maximal radius at which the parameterization in radial space — the 
Fourier transform of the parameterization in momentum space — allows for a somewhat 
free variation of the density. For the nuclei of interest here AQ ~ 0.2 fm“^. While the 
numerical values of the extracted form factors no longer show the statistical fluctuations 
of the data, the error bars do account for both the random and the systematic errors 
of the data. These error bars are much more quantitative than the fluctuations of the 
usual form factors, which in general are taken as a visual “measure” of the accuracy of 
the form factors. 

In order to determine the form factors we have htted the world data using the highly 
flexible sum-of-gaussians (SOG) parameterization [13]. The main restriction, introduced 
by the SOG in r-space, is the rms-radius of the gaussians, which is chosen to be well 
below the rms-radius of the proton. The resulting form factors and their error bars are 
shown in the hgures of section |H 

The approach described above may seem ’unconventional’, but it is perfectly 
analogous to the procedure employed since decades in nucleon-nucleon (NN) scattering. 
The NN data (cross sections and analyzing powers) are no longer used directly, but the 
information is condensed in the phase-shifts extracted from a global ht to the data. 
The procedure we employ corresponds to the standard energy-dependent phase-shift 
analyses m. There are however two main differences: (i) we do not prune the data, 
i.e. eliminate some 30% of the data in order to get the y^/datum down to ~ 1 mi; (ii) 
we do take into account the absolute normalization of the data, rather than floating the 
cross sections in order to produce a y^/datum~ 1. 

The hts to the world data also provide the most accurate values for integral moments 
such as rms-radii and the Zemach moments [IHl [16] of interest for the interpretation of 
transition energies in atomic nuclei. The radii are also quoted in section 01 The moments 
in general have been obtained by constraining the shape of the large-r density to the one 
expected from our understanding of nuclear wave functions; at large r they must fall 
like a Whittaker function depending on the nucleon separation energy. This constraint 
is needed as the determination of e.g. the rms-radius involves an extrapolation from the 
Q-region sensitive to hnite size — typically 0.5 fm' Q <1.2 im. ^ — to Q = 0 where 
the radii are extracted. The difficulties of this (implicit) extrapolation are discussed in 
Ref. HI]. 
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2.2. Two-photon corrections 

During the last years, it has become clear that two-photon exchange beyond the one 
accounted for via the Coulomb corrections (exchange of a soft photon in addition to hard 
photon responsible for the scattering) can contribute. At large momentum transfer the 
exchange of two hard photons can become important, as was shown in particular for 
the proton, where the two-photon corrections resolved the discrepancy between L/T- 
separation and polarization transfer results for the proton charge form factor. For a 
review see Ref. [I8] . 

Two-photon exchange is difficult to calculate due to the complexities introduced by 
the intermediary states, which can be the nucleus of interest in its ground state or in any 
excited state. Intermediary excitation of individual nucleons can also contribute. For the 
deuteron, calculations are available by Kobushkin et al. [I9] and Dong et al. [2011211122]. 
While Dong et al. calculated the contribution where the exchange takes place with only 
one nucleon, Kobushkin et al. assume that the exchange involves both nucleons. For 
^H, a calculation is available from Kobushkin and Timoshenko [23]. For the deuteron, 
the calculated two-photon exchange contributions are of the order of 10 % of the typical 
differences between theory and experiment, for A = 3 they can, at selected places, 
amount to up to 30 %. Since the contributions of two-photon exchange at the present 
time are still quite uncertain, we will not include them. 

2.3. Experiments on deuteron 

A large number of experiments on elastic electron-deuteron scattering have been 
performed since ~1960, producing a total of some 500 data points [ 21 ] - [ 19 ]. We 
discuss below only a selection which was particularly important in hxing the deuteron 
form factors. For the analysis of the data in terms of form factors all data will be 
included. 

The most accurate data at low momentum transfers come from the experiments of 
Simon et al. m performed at Mainz, in the range Q = 0.2 — 2 fm the experiment 
of Platchkov et al. [H] performed at the Saclay ALS, with Q = 1.2 — 4.2 fm“^, and 
the experiment of Berard et al. EH carried out at the Monterey accelerator, with 
Q = 0.2 —0.7 fm~^. These experiments reached accuracies of the order of 1%, employing 
liquid deuterium targets and high-resolution magnetic spectrometers for the electron 
detection. In some special cases gas targets of better known thickness and hxed-angle 
spectrometers of well controlled solid angle were used as a supplement to achieve better 
accuracy on absolute cross sections. 

Data at medium momentum transfers come, among others, from two experiments 
carried out at JLab by Alexa et al. [25], for Q = 4.2 — 12.2 fm“^, and Abbott et al. [2l] . 
for Q = 4 — 6.6 fm“^. Although these experiments quote rather small errors, the data are 
not entirely consistent; we suspect that, in the early days of JLab operation, the beam 
energies were not very accurately known. A more recent JLab experiment of Bosted et 
al. [50], in the range Q = 3.7 — 5.7 fm~^, could not resolve the difference in favor of one 
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of the experiments. 

Data up to the highest momentum transfers were also measured at SLAG by Arnold 
et al. IZ71, in the range Q = 4.5 — 10.1 fm This experiment, as the ones performed 
at JLab, detected the scattered electron and recoil deuteron in coincidence, such as to 
cleanly identify elastic scattering despite insufficient energy resolution. 

An important class of experiments provided the backward-angle data needed to 
determine the magnetic form factor. Although these experiments provided less accurate 
cross sections, they are most valuable because they are generally totally dominated 
by Gm{Q) and thus less dependent on the (error-enhancing) L/T-separation. The 
experiments of Benaksas et al. ED], with Q = 1.7 — 2.2 fm and Ganichot et al. mi, 
with Q = 0.7 — 2.4 fm”^, were carried out at 180° scattering angle at the Orsay 
accelerator, the data of Rand et al. j46], with Q = 2.2 — 3.1 fm“^, come from a 180° 
experiment done at the Stanford HEPL machine. Precise large-angle data were also 
provided by the experiment of Auffret et al. [29], with Q = 2.4 — 4.2 fm“^, with a 
scattering angle of 155° at the Saclay ALS. 

A special class of data, only recently accessible to experiment, involves measurement 
of tensor polarization observables, accessible with deuterons with spin aligned in the 
direction of the momentum transfer. When working with a tensor-polarized deuteron 
target, or when measuring the tensor polarization of the recoiling deuteron, the quantity 
T2o{Q) can be obtained (the other tensor observables T 2 i{Q) and T 22 {Q) are not that 
useful). With the knowledge of T 2 o{Q), which basically depends on the Gq{Q)/Gc{Q) 
ratio [see equation (II.bh ]. it becomes possible to separate the form factors Gc{Q) and 
Gq{Q) which cannot be separated via cross-section measurements alone. In the Q-range 
of interest here, Gc{Q) is of particular interest due to the presence of a diffraction 
feature, which is very sensitive to the ingredients of the theoretical calculations. 

Measurements of T 2 o(Q) have been performed at various laboratories [5T] - [60] . 
They have been performed both at storage rings using internal, polarized deuteron 
atomic beams, or with external electron beams using polarimeters to detect the recoil 
tensor polarization. The most extensive set of data, measured by Abbott et al. at 
JLab [5^, comes from an experiment that used a polarimeter based on the analyzing 
power of the p{d, 2p)n reaction which had been calibrated using a polarized deuteron 
beam from the Saturne accelerator. Today, data are available up to momentum transfers 
of 6.6 fm“^. 

For the region Q < 7 fm“^, where the Gc/f^M-separation can be performed, the 
data base comprises some 492 data points. The ht yields a of 549 when taking into 
account only the statistical errors of the data. 

To make a connection between the present global analysis and the previous 
measurements discussed above, we show in £gure[T]the ratio of the experimental values 
for the usual structure functions A{Q), B{Q), T 2 o{Q) existing in the literature to the 
ones obtained in the present global analysis. Note, however, that in our global analysis 
we use roughly twice as many data points, since there are many data available that have 
not led to a determination of the deuteron structure functions. It also should be noted 
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Figure 1. (Color online) Ratio of the structure functions A{Q), B{Q), T 2 o{Q) as 
measured by the different experiments to the presented global analysis (GA). The 
different labels to identify the experiments correspond in order, for A{Q), to Refs. [Ml 
|30l[32l[39l[22lll7l|45l[25l[24], for B{Q) to Refs. |38l [Ml EHl SSI Eg li Ell ESI ESI [28] , 
and for T 2 q{Q) to Refs. miESlESlEilESlESlESlEelEllEOlEE]. Note that Refs. [54] 
and [59] refer both to the label “Bates (1994)”. 
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that the structure functions from the literature have been derived without accounting 
for Coulomb distortion. 

2.4- Experiments for 

The nucleus represents a particular challenge due to its radioactive nature (with a 
half-life of 12 years); experiments involve up to 10^ Curies of material. This requires 
strict safety measures, particularly since high-intensity beams as available at modern 
accelerators can easily melt a hole into the windows of the target container. 

The earliest experiment was performed at the Stanford HEPL accelerator by Collard 
et al. [6T] using a high-pressure (100 bar) gas-target. This experiment reached Q = 2.8 
fm“^. Another experiment with a gas target was performed about 20 years later at the 
MIT/Bates accelerator. This experiment, by Beck et al. [62], used tritium stored in 
an Uranium-oven; in the target the tritium was cooled down to 45 K, thus allowing to 
run at a more modest pressure of 15 bar. This experiment reached a similar maximal 
momentum transfer as the previous one, but better statistical accuracy due to the 
achieved higher luminosity. 

The most extensive set of data comes from the experiment performed at the Saclay 
ALS accelerator by Amroun et al. [6^ l64] . In order to reach a much higher luminosity, 
this experiment used liquid tritium, cooled to a temperature of 20 K. With a novel 
target system — cooled 5 cm long cylindrical target plus warm storage vessel of 0.2 liter 
volume — a permanently sealed system could be employed with most of the tritium 
(98%) in the target at a pressure as low as 3 bar when the target was in operation. This 
led to the highest luminosity with only 10"^ Curies of tritium. This experiment reached a 
maximum momentum transfer of Q = 5 fm“^ and covered the region of the diffraction 
minimum and maximum expected from the then already known data on ^He. 

All three experiments provided data at both forward and backward angles, 
collecting about 190 data points, so that a Rosenbluth separation into charge and 
magnetic form factors is possible. This separation has been performed as described 
above. The of the SOG-ht is 340, mainly due to a difference in normalization of 
the Saclay/Bates data sets. 

2.5. Experiments on ^He 

A good target for ^He represents a lesser challenge than an ^H-target, although still a 
major effort is needed to reach an adequate target thickness and luminosity. 

Data at low momentum transfers have been measured by Szalata et al. [65] using 
the NBS accelerator {Q < 0.6 fm“^), vonGunten [66] at the Darmstadt machine, and 
Ottermann et al. EH at Mainz {Q < 1.9 fm ^). The experiment of Dunn et al. [68] . 
performed at the Bates laboratory, provided data up to Q ~3.3 fm“^. 

The region of the predicted diffraction minimum and maximum of the charge form 
factor was reached by the experiment of McCarthy et al. [69] [TO] . Contrary to all other 
experiments, McCarthy et al. used a liquid ^He-target, cooled by superfluid ^He. This 


Electromagnetic Structure of Few-Nucleon Ground States 


12 



1.6 



1 *^ - 

xl 


+ Saclay(1994) 

X Saclay{1982) 

+ Bates (1983) 

O Bates (1987) 

• NBS(1977) 

□ Mainz (1985) 

A HEPL (1970- 1977) 

^ SLAC(1978) 

▼ HEPL (1965) 

► Bates (2001) 

I ^ • Darmstadt (1982) 

6 7 


Q [fm"'] 


Figure 2. (Color online) Ratio of experimental cross section over the one obtained 
by the present global analysis for ^He. The different labels to identify the experiments 
correspond in order to Refs. [SainilMlIMlIMlETlIMlIIQlEaielllzillM]- 


produced a much higher target thickness and luminosity, and allowed to measure very 
small cross sections, leading to a maximum momentum transfer of 4.5 fm“^. Use of 
liquid ^He came at the expense of a rather complicated target setup. 

The diffraction feature in the magnetic form factor was observed in the experiment 
of Cavedon et al. 1711 performed at the Saclay ALS accelerator. Additional backward- 
angle points were measured at 160° at Bates by Nakagawa et al. [73]. With this 
experiment, the data on the magnetic form factor now reach a maximum transfer of 
6.5 fm“h 

Data up to the highest momentum transfers, Q ~10 fm“^, were provided by Arnold 
et al. na with an experiment carried out at SLAG at forward scattering angles. With 
a high-pressure gas target cooled by liquid hydrogen, this experiment could reach the 
extremely small form factors near the second diffraction minimum. Due to the limited 
energy resolution of the spectrometers, scattered electron and recoil ^He had to be 
detected in coincidence in order to cleanly identify elastic scattering. 

The data set for ^He comprises some 310 data points; the ht has a of 347 when 
ignoring the systematic errors of the data. The ratio of experimental cross section over 
the one obtained by the present global analysis for Q < 7 fm“^ is shown in hgure O 
Note that an experiment going to similarly large Q values has been performed at JLab, 
but the data have not yet been published. 

For the A = 3 system, it is of interest to study, besides the form factors of the 
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individual nuclei, the isoscalar and isovector combinations as well. In particular, the 
isoscalar charge form factor, as defined in equation fll.26p . is more directly comparable to 
the charge form factors of deuteron and “^He. These isoscalar and isovector combinations 
can be easily calculated, together with their uncertainties, from the form factors fitted 
to the experimental data. 

2.6. Experiments on '^He 

This nucleus is particularly tightly bound, resulting in the largest density in the inner 
region. This property makes the form factors of '^He especially interesting. Due to its 
Jj = 0 nature, the determination of the form factor is easier than for all other light 
nuclei, and a total of 190 data points are now available. 

Experiments at low momentum transfers have been carried out by vonGunten [66] 
and Erich et al. m at Darmstadt, and by Ottermann et al. [67] at Mainz, the 
latter reaching Q ~2 fm~^. The data at medium momentum transfers come from 
two experiments performed at the Stanford HEPL machine by Frosch et al. m and 
McCarthy et al. 1701. These experiments, performed using liquid targets, covered the 
region up to Q =4.5 fm“^ and thus the region of the first diffraction minimum and 
maximum. 

The highest-Q data were measured by Arnold et al. 1721 at SLAG, reaching Q = 7.9 
fm“^. The recent publication by Camsonne et al. [76] provided data up to Q = 8.8 
fm“^. Profiting from the high luminosity achievable at JLab, this experiment covered 
the region of the second diffraction feature. In the region of overlap with the data of 
Ref. [72] some disagreement is visible, whose origin, at the present time, has not been 
cleared up. 

3. Theoretical approaches 

In this section we discuss the different theoretical frameworks adopted to study A = 2—4 
electromagnetic form factors. Section [3T] reviews the conventional approach, which uses 
phenomenological realistic models for the nuclear interactions and currents, section 13.21 
reviews the chiral effective field theory approach, in which both currents and interactions 
are consistently derived from chiral Lagrangians, and section 13.31 reviews results from 
the relativistic covariant spectator theory, where (using the language of covariant field 
theory) currents consistent with the relativistic dynamics are constructed. Finally, 
section 13.41 presents an appraisal and comparison of these various methods. 

3.1. The conventional approach 

The conventional approach views the nucleus as made up of nucleons interacting among 
themselves via two- and many-body potentials, and with external electroweak fields 
via one- and many-body currents. It assumes that all other sub-nucleonic degrees 
of freedom, involving, for example, the excitation of nucleon resonances such as A 
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isobars, can be eliminated in favor of these effective potentials and currents, acting 
on nucleon coordinates. (For the size of effects of single A-isobar excitations in two- 
and three-nucleon form factors see, for instance, the pioneering work by Sauer and 
collaborators m -|80].i The validity of this greatly simplified description, in which 
color-carrying quarks and gluons (the degrees of freedom of quantum chromodynamics, 
QCD) are assembled into colorless clusters (the nucleons), and these clusters are taken 
as effective constituents of the nucleus, ultimately rests on the success it has achieved 
in the quantitative description of many nuclear properties. 

3.1.1. Nuclear potentials In the current version of the conventional approach, the 
nuclear Hamiltonian is taken to consist of a non-relativistic kinetic energy term and 
of two- and three-body potential energy terms. 



(3.1) 


i i<j i<j<k 


The kinetic energy term is predominantly charge-independent, though it has a small 
charge-symmetry-breaking (CSB) component due to the difference between proton and 
neutron masses which leads to a very small correction to nuclear energies [1]. The two- 
body potential consists of a long-range part, for inter-nucleon separation r >2 fm, due 
to one-pion exchange (OPE), and intermediate- and short-range parts, for, respectively, 
1 fin < r < 2 fm and r < 1 fm, which are derived from theoretical arguments and 
are constrained by fits to nucleon-nucleon (NN) elastic scattering data typically up to 
lab kinetic energies of 350 MeV, slightly above the pion production threshold. These 
potentials are customarily referred to as “realistic potentials” in the literature. In the 
case of the CD-Bonn (CDB) momentum-space (and strongly non-local) potential [8T] . 
the short- and intermediate-range parts are parametrized in terms of p and cu vector- 
meson exchanges as well as the exchange of two effective scalar mesons, whose masses 
are about 450 (350) MeV and 1220 (800) MeV, as determined by hts to pp data {np 
data in isospin channel T = 0). In the case of the Argonne Uig (AV18) conhguration- 
space potential [82], the intermediate range part is parametrized in terms of two-pion 
exchange (TPE), based on, but not consistently derived from, a held-theory analysis 
of (direct and crossed) box diagrams with intermediate nucleons and A isobars [83] . 
while its short-range part is represented by spin-isospin (and momentum-dependent) 
operators multiplied by Woods-Saxon radial functions. This potential is used in most 
of the results obtained in the conventional approach and presented in later sections of 
this review; for details see Ref. [82] . 

Modern realistic potentials contain isospin-breaking (IB) terms. At the level of 
accuracy required, electromagnetic interactions, along with strong interactions, have to 
be specified in order to ht the data precisely. These electromagnetic interactions consist 
of one- and two-photon Coulomb terms as well as Darwin-Foldy, vacuum polarization, 
and magnetic moment contributions [H]. The full potential Vij is then the sum of isospin- 
conserving strong-interaction terms vj^, specified electromagnetic-interaction terms nh 
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up to order where a is the fine structure constant, and finally additional isospin¬ 
breaking strong-interaction terms 

In the context of the conventional framework described here, it is an established 
fact that accurate calculations based on realistic two-nucleon potentials find that the 
observed energy spectra of light nuclei with mass number A in the range 3 < A < 12 
are under-predicted [85l [86] and that the empirical saturation density of nuclear matter 
is over-predicted by roughly a factor of 2 [871 EH]- In the specific case of the few- 
nucleon systems of interest in this review, (essentially exact) hyperspherical-harmonics 
calculations lead to triton binding energies ranging from about 7.6 MeV for a local 
potential like the AV18 to about 8.0 MeV for a strongly non-local potential like the CDB, 
compared with the experimental value of 8.48 MeV. For ^He, the binding energy ranges 
from 24.0 MeV for the AV18 to 26.3 MeV for the CDB, compared to the experimental 
value of 28.3 MeV [85] . 

Several effects could be important in reproducing the binding energies of nuclei. 
Two of them are immediately apparent: relativistic corrections and three-nucleon 
interactions. It has long been known that these effects cannot be completely separated, 
i.e. that they are theoretically related (for a brief discussion see Section lT4l below). 
Furthermore, their contributions are comparable. In the “^He ground state, for example 
the (non-relativistic) kinetic energy is about 100 MeV, and one would expect 1-2 % 
relativistic corrections of this value, amounting to 1-2 MeV. Three-nucleon potentials 
lead to corrections of similar size. At long range, the three-nucleon potential is of the 
well known Fujita-Miyazawa type [89], corresponding to single-pion exchanges between 
three nucleons with intermediate excitation of a A-isobar resonance. In coupled-channel 
calculations of the three-nucleon bound state in which the A-isobar is treated as an 
active degree of freedom, this three-nucleon force yields almost 1 MeV of additional 
binding energy, but it is to a large part cancelled by the dispersive contribution of 
intermediate A-excitations to two-nucleon scattering [90l|911[78]. In ^He, the presence of 
this relatively low-lying resonance produces a three-nucleon potential whose contribution 
is of roughly a few MeV. 

The Fujita-Miyazawa three-nucleon potential has the following structure 


V.% = 


1 {Tj ■ T? ’ Ti ■ '^k\ T ^ 1 [Ti ' D ’ Ft ' 


(3.2) 


where Xij = Y.,,(rij) cTi ■ cTj + T.„{rij) Sij, {Xij,Xik} {[Xij,Xik]) is the anticommutator 
(commutator), and the YTr{r) and Tj^ir) radial functions are those entering the OPE two- 
nucleon potential. In the approach developed in Ref. [92] and adopted in the calculations 
of few-nucleon form factors reported in this review, the attractive three-nucleon potential 
above is supplemented by a purely central short-range term, i.e., 

V.it = + V,% . V,% = U„Y, . (3,3) 


The term is of two-pion-exchange range on each of the two legs. It is meant 
to simulate the dispersive effects that are required when integrating out A degrees 
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of freedom. These terms are repulsive and are taken to be independent of spin and 
isospin. The constant ^ 2 ^ in and Uq in in combination with the AV18 two- 
nucleon potential, are adjusted to reproduce the triton binding energy and to provide 
additional repulsion in hypernetted-chain variational calculations of nuclear matter near 
the equilibrium density. The three-nucleon potential in equation (13.31) is denoted as the 
Urbana IX (UIX) model in what is to follow, and AV18/UIX is used to denote the 
Hamiltonian including the AV18 two-nucleon and UIX three-nucleon potentials. 

Before moving on to a discussion of the nuclear electromagnetic current operator, it 
should be pointed out that recent developments in quantum Monte Carlo methods and, 
in particular. Green’s function Monte Carlo methods have made it possible to carry out 
essentially exact calculations of the energy spectra of low-lying states of light (s- and 
p-shell) nuclei in the mass range up to A = 12 (^^C). These calculations have exposed 
the inadequacy of the AV18/UIX model to satisfactorily reproduce the observed spectra 
of A=6-12 nuclei, and have led to the development of a new model for the three-nucleon 
potential—the Illinois-7 (IL7) model [93]. The latter incorporates V^l, but also includes 
terms involving multi-pion exchanges and intermediate A’s as well as a representation of 
the short-range term which now retains isospin dependence. It is characterized by 
three parameters, which have been determined, in combination with the AV18, by htting 
the low-lying states of A=3-10 nuclei. The resulting AV18/IL7 Hamiltonian then leads 
to predictions of several ground- and excited-state energies, including the p—^He elastic 
scattering observables [9l| and the ^^C ground- and Hoyle-state energies [93], in very 
good agreement with the corresponding empirical values. However, in the few-nucleon 
systems of interest here, bound-state energies obtained with either the AV18/IL7 or 
AV18/UIX Hamiltonian are not signihcantly different. 


3.1.2. Nuclear electromagnetic charge and current operators A fundamental aspect in 
the description of electromagnetic (and weak) processes in nuclei is the derivation of 
a consistent set of nuclear electromagnetic (and weak) currents. The leading terms 
are expected to be those associated with the charges and convection currents of the 
individual protons, and the spin-magnetization currents of the individual protons and 
neutrons, which in conhguration space are 


Pi(q) = 


^/Y^^~Qy4rn3 


ei(Q) - ^ [2 EiiQ) - g{Q)] q ■ (ct* X Pi) 




(3.4) 


ji(q) = ■ (3.5) 

2m 2m 

and follow from a non-relativistic expansion of the covariant single-nucleon current, 
including corrections to the operators up to order jvn?. Here q and cu are as defined 
in section [H is the momentum operator of nucleon i with its charge and magnetization 
distributions described by the form factors Ciiff) and 


-i{Q) = 2 [G|W) + GpQ) r,J , 


(3.6) 
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^ [GmW) + G!:,W) r, J . (3,7) 

where G^{Q) and Gf^{Q), and G^{Q) and G\^{Q), are, respectively, the isoscalar 
electric and magnetic, and isovector electric and magnetic, combinations of the proton 
and neutron form factors, normalized as Gf(0) = G^(0) = 1, Gff(O) = /i'^, and 
G\j{fS) = /i'^, with and jjX denoting the isoscalar and isovector combinations of the 
proton and neutron magnetic moments, = 0.88 and jjX = 4.706 in units of nuclear 
magnetons. These form factors are obtained from hts to elastic electron scattering data 
off the proton and deuteron [95] . 

The current and charge operators given in equations (13.dh and (13.51) are usually 
referred to as impulse approximation (lA). There is ample evidence that these lA 
operators are inadequate, especially for the description of isovector currents. This 
evidence comes from studies of a variety of photo- and electron-nuclear observables 
at low and intermediate values of energy and momentum transfers, especially in light 
nuclei {A < 12) for which essentially exact calculations can be carried out. Experimental 
data are poorly reproduced in lA. Well known examples are the, classic by now, 10% 
underestimate of the n — p radiative capture cross section at thermal neutron energies, 
which in fact provided the initial impetus to consider two-body terms in the nuclear 
current operator [96]- [99], the 15% underestimate of the isovector magnetic moment of 
the trinucleons, the large discrepancies between the experimental and calculated charge 
and magnetic form factors of the hydrogen and helium isotopes, particularly in the 
hrst diffraction region at momentum transfers in the range of 3-3.5 fm“^, the large 
underestimate of the n — d and n—^He radiative captures [1001 llOlj . and, hnally, the 
signihcant underestimate, in some cases even of 40%, of magnetic moments and Mi 
radiative transition rates in A=7-9 nuclei una. 

Many-body terms in the nuclear electromagnetic charge and current operators 
arise quite naturally in a meson-exchange picture or when the excitations of nucleon 
resonances, such as the A isobar, are taken into account. There is a very large body 
of work dealing with the problem of constructing these electromagnetic many-body 
operators from meson-exchange theory, and we defer to a number of reviews [10311104] 
for a summary of efforts along those lines. Here we will describe an approach, originally 
proposed by Riska [96l [971 EH] , that leads to conserved currents, even in the presence 
of two- and three-nucleon potentials, not necessarily derived from meson-exchange 
mechanisms (as is the case for the AV18 and UIX models). This approach has been 
consistently used by the ANL/JLab/LANL/Pisa group to study many electromagnetic 
processes in light nuclei (up to ^^C), and has proved to be quite successful in providing 
predictions systematically in close agreement with experiment. 

The dominant part of any realistic NN potential consists of static (momentum- 
independent) terms, and leading electromagnetic two-body charge and current operators 
are derived from these terms, specihcally the isospin-dependent central, spin, and tensor 
components. The latter are assumed to be due to exchanges of effective pseudo-scalar 
(PS or TT-like) and vector (V or p-like) mesons, and the corresponding charge and current 
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operators are constructed from non-relativistic reductions of Feynman amplitudes with 
the TT-like and p-like effective propagators. For the vr-like case (we defer to Ref. |6] for 
a complete listing), they read 


j.y(k.,k,) =iGUQ){ 


V ( 


Ti X T, 


3)^ 


Vps{kj 


kj — k. 


■ ki 




(k„ k,) = [Ff (g) T. ■ T, + Fr (Q) r.,] 


vpsjkj 
2 m 


(Ti ■ q CT ■ k, + i 


3 


(3.8) 


(3.9) 


where kj and kj are the fractional momenta delivered to nucleons i and j, with 
q = kj + kj, and Vps{k) is projected out of the (isospin-dependent) spin and tensor 
components of the potential (see equation (2.26) of Ref. |100] ). The Dirac nucleon 
electromagnetic form factors F ^' are related to those introduced previously via 
+ rj /(1+p) with p = Q ‘^/(4 m^), and therefore differ from by 

relativistic corrections proportional to rj. By construction, the longitudinal components 
of the resulting and jh currents satisfy current conservation with the static part of 
the AV18, denoted Ujj (static) below. 


q • [jps(q) + jn(q)] = [Di(static) , Pi(q) + Pj{q) ] . (3.10) 

Hence the use in equation fl3.8l) of the form factor G'^{Q) entering Pi(q). Of course, 
the continuity equation poses no restrictions on transverse components of the current, 
in particular on electromagnetic hadronic form factors that may be used in these 
components. Ignoring this ambiguity, the choice is made here for both longitudinal 
and transverse components. 

Additional conserved currents follow from minimal substitution in the momentum 
dependent part of Vij, denoted as (non-static). This momentum dependence enters 
explicitly via the spin-orbit, quadratic orbital angular momentum, and quadratic spin- 
orbit operators, and implicitly via Tj • Tj, which for two nucleons can be expressed in 
terms of space- and spin-exchange operators as in 

Tj • Tj = -I - {I + (Ti ■ aj) . (3.11) 


Both the explicit and implicit (via Ti-Tj) momentum-dependent terms need to be gauged 
in order to construct exactly conserved currents with Ujj(non-static). The procedure, 
including the ambiguities inherent in its implementation, is described in Ref. |100j . In 
contrast to the purely isovector and jh, the currents from Ujj(non-static) have both 
isoscalar and isovector terms, which, however, due to their short-range nature, lead to 
contributions that are typically much smaller (in magnitude) than those generated by 
and jh. 

Finally, conserved three-body currents associated with the term of the Vijk 
have also been derived by assuming that this term originates from the exchange of 
effective PS and V mesons with excitation of an intermediate A isobar Hffll. However, 
their contributions have been found to be generally negligible, except for some of the 
polarization observables, like T 20 and T 21 , measured in the proton-deuteron radiative 
capture at low energy HOT]. 
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It is important to stress that the two- and three-body charge and current operators 
discussed so far have no free parameters, and that their short-range behavior is consistent 
with that of the potential, which is ultimately constrained by NN scattering data. 
It is also worthwhile noting that the two-body charge operators vanish at vanishing 
momentum transfer, as they must in order to conserve the overall charge of the nucleus. 
They also vanish in the static limit (m —)■ cxo). It was pointed out by Friar |105] long 
ago that a proper derivation of the leading two-body charge operator of pion range, the 

in equation (13.Qh . necessarily entails the study of non-static corrections to the OPE 
potential, and that in particular its form depends on the specihc, but arbitrary, off- 
the-energy shell extension—that is, on the corrections beyond the static limit, such as 
those induced by retardation effects—adopted for it. Furthermore, he showed that these 
different operators and corresponding (non-static) OPE potentials are related to each 
other by a unitary transformation, which implies that their intrinsic lack of uniqueness 
has no consequence for physical observables. However, non-static corrections to the OPE 
term are not considered in the AV18. It is also reassuring to note that a pion-range 
charge operator of the form given in equation fl3.9l) has been derived in the context of 
chiral effective held theory; see below. 

Additional short-range isoscalar (isovector) two-body charge and current operators 
follow from the pvry (uTrj) transition mechanism. The pvry and UTrj currents are purely 
transverse and therefore unconstrained by current conservation. The coupling constants, 
and hadronic and electromagnetic form factors at the pNN, uNN, pvry, and cuvry 
vertices are poorly known [6]. However, with the exception of the pvry current which 
gives a signihcant contribution to the deuteron magnetic form factor, generally these 
operators lead to very small corrections to the charge and magnetic form factors of the 
few-nucleon systems of interest in this review. 

Finally, there are purely transverse many-body currents arising from Mi-excitation 
of A resonances. They have been derived in a number of different approaches, the most 
sophisticated of which is based on the explicit inclusion of A isobar degrees of freedom in 
nuclear wave functions. In this approach, known as the transition-correlation-operator 
(TCO) method and originally developed in Ref. |106] . the nuclear wave function is 
written as 


4'k+a= « 


(3.12) 


. i<j 


where T is the purely nucleonic component and S is the symmetrizer. The transition 
operators convert NN into NA and A A pairs and are obtained from two- 

body bound and low-energy scattering solutions of the full N + A coupled-channel 
problem |in6] . including transition potentials vfj^{NN NA) and v]j^{NN —)■ AA). 
Indeed, the simpler perturbative treatment of A-isobar degrees of freedom, commonly 
used in estimating the A-current contributions, uses the approximation 



m — niA 


1 


[u5^(AA ^ NA) + i-j] + 
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and mA (1232 MeV) is the A mass. This perturbative treatment has been found 
to be inappropriate, since it overestimates A contributions. In the presence of an 
electromagnetic held, N ^ A and A —)■ A couplings need to be accounted for, and these 
couplings (and associated electromagnetic form factors) are taken from N{e, e') data in 
the resonance region. In practice, the currents arising from A resonance excitation can 
be reduced to effective two- and many-body operators depending on f/™, but acting 
only on the nucleonic component T of the full wave function. The TCO method is used 
in some of the calculations reported in the present review. 

In the conventional framework based on instant-form Hamiltonian dynamics, it is 
possible to perform calculations within a n/c expansion scheme, in which the Poincare 
covariance of the theory is satished to order (n/c)^ |1071110811109] . In this approach the 
many-body Hamiltonian is written as 

^ = E(\/p? + -™) + E+ E [Vi,t + SVtMPm)] . (3.14) 

i ' ^ i<j i<j<k 


where the relativistic expression for the kinetic energy is used, and Vij and Vijk are, 
respectively, the two- and three-nucleon potentials in the corresponding rest frames, 
while the so-called boost corrections 6vij(Pij) and SVijkiPijk) depend on the total 
momenta Pij and Pijk of, respectively, the two- and three-body subsystems, and vanish 
when Pij = 0 and Pijk = 0. These boost corrections are related to the rest-frame 
potentials by the requirement that the commutation relations of the Poincare group be 
satished to order (n/c)^ |1091 IllOj . The effects of these boost corrections on the binding 
energies of light nuclei have been studied in Refs. |1101 l 86 ] . In these calculations, 
the rest-frame two-body potential Vij must be rehtted to NN data, and obviously the 
parameters present in Vijk must also be recalibrated. The results of a comparison with 
a phase-equivalent nonrelativistic Hamiltonian show that the relativistic corrections to 
the binding energies are repulsive: in ^He, for example, they amount to about 2 MeV. 

The approach above has been used to carry out a deuteron form factor 
calculation m, by consistently including also the (n/c)^ corrections arising from the 
boosting of the deuteron wave function. 


V'v(p) = ^ 
^/7 L 


4m 


[CTi - 0-2) X P 


V'o(P||/7,P±) 


(3.15) 


where v is the velocity of the moving frame, 7 = I/a/I — p|| and p_L are the 

components of the momentum p parallel and perpendicular to the velocity v, and fjo 
is the rest frame wave function. The calculation of Ref. m also retained the full 
covariant structure of the one-body current and two-body pvry current, and the pion- 
like and p-like two-body terms (correct to order (n/c)^) that contribute to the isoscalar 
charge operator. Some of the results of this approach are presented below. 


3.2. The chiral effective field theory approach 

The last two decades have witnessed signihcant developments in nuclear chiral effective 
held theory (yEFT), originally proposed by Weinberg in a series of papers in the early 
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nineties |112l I113[ I114j . The (approximate) chiral symmetry exhibited by quantum 
chromodynamics (QCD) severely restricts the form of the interactions of pions among 
themselves and with other particles. In particular, the pion couples to baryons, such as 
nucleons and A-isobars, by powers of its momentum Q, and the Lagrangian describing 
these interactions can be expanded in powers of Q/A^, where ~ 1 GeV specifies 
the chiral-symmetry breaking scale. In what follows, Q will represent not only the 
momentum of the pion, but may also the generic value of the momentum of other 
particles. It is assumed to be less or equal to the pion mass. As a result, classes of 
Lagrangians emerge, each characterized by a given power of Q/A^ and each involving 
a certain number of unknown coefficients, the so called low-energy constants (LEC’s). 
These LEC’s are then determined by fits to experimental data (see, for example, the 
review papers pg and pg, and references therein). Thus, yEFT provides, on the 
one hand, a direct connection between QCD and its symmetries, in particular chiral 
symmetry, and the strong and electroweak interactions in nuclei, and, on the other 
hand, a practical calculational scheme, which can, at least in principle, be improved 
systematically. In this sense, it can be justifiably argued to have put low-energy few- 
nucleon physics on a more fundamental basis. 

Within the nuclear xEFT approach, a variety of studies have been carried out 
in the strong-interaction sector dealing with the derivation of two- and three-nucleon 
potentials [117j - [125] and accompanying isospin-symmetry-breaking corrections [126] 
- ra, and in the electroweak sector dealing with the derivation of parity-violating 
two-nucleon potentials induced by hadronic weak interactions [130] - [133] and the 
construction of nuclear electroweak currents [134] . In this review, the focus is on nuclear 
electromagnetic charge and current operators. These were originally derived up to one 
loop level in the heavy-baryon formulation of covariant perturbation theory by Park et 
al. ina. More recently, however, two independent derivations, based on time-ordered 
perturbation theory (TOPT), have appeared in the literature, one by some of the present 
authors [1361113711138] and the other by Kolling et al. [13911140] . In the following, we 
outline the derivation of these operators, deferring the discussion of some of the more 
technical aspects to the original papers [13611137] 1138] . 

3.2.1. Interaction Hamiltonians In the simplest implementation, yEFT Lagrangians 
are constructed in terms of nucleon and pion degrees of freedom. By now, the procedure 
by which this is accomplished has been codihed in a number of papers m, and 
ttN and tttt Lagrangians, denoted respectively as and C^\ have been derived 
up to high order in the chiral expansion. Contributions arising from the inclusion of 
additional degrees of freedom, such as A-resonances and heavier mesons, are effectively 
subsumed in the LEC’s entering these and Lagrangians. In principle, they 
contain an inhnite number of interactions compatible with the QCD symmetries, but 
as mentioned above, the transition amplitudes obtained from them can be expanded 
in powers of Q/A^, and at each given order of the expansion, the number of terms 
contributing to the amplitude is hnite |112[ 1113] 1114] . The canonical formalism is 




















































Electromagnetic Structure of Few-Nucleon Ground States 


22 


used to construct the Hamiltonians from the chiral Lagrangians. Only terms entering 
the two-nucleon potential and electromagnetic charge and current operators up to one 
loop are listed below. As it will become apparent later, this requires the potential up to 
order (Q/A^)^ and the electromagnetic operators up to order (Q/A^)^ for a two-nucleon 
system. However, it is worthwhile pointing out that at the present time two-nucleon 
potentials have been derived and widely used up to order (Q/A^)^ (requiring two-loop 
contributions). Very recently, a new derivation up to order (Q/A^)® has appeared |142] . 
Indeed, some of these high-order potentials have been used, in conjunction with the 
one-loop electromagnetic operators, in the calculations of static properties and form 
factors of A=2 and 3 nuclei reported below. From a chiral counting perspective, these 
calculations are not consistent, since strict adherence to consistency would require going 
up to order (Q/A^)^ in the derivation of the electromagnetic operators, a rather daunting 
task. It is also unclear at this point how many new LEC’s would enter, in addition to 
the hve at order (Q/A^)^ (see below); if there were to be too many of them, this would 
obviously reduce substantially the predictive power of the theory, since there are only 
a limited number of electromagnetic observables in the few-nucleon systems (including 
single nucleons) to constrain these LEC’s. 

Setting aside these considerations, the subset of interaction terms in and 

in the ttN sector, and in the tttt sector, relevant to the derivation of potentials 
and electromagnetic operators at one loop level leads to the following Hamiltonians 


H. 
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= /dxVt 
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N , (3.19) 


where qa, fw, e, and m are, respectively, the nucleon axial coupling constant, pion 
decay amplitude, proton electric charge, and nucleon mass, and the parameters d[ are 
(yet to be determined) LEC’s. The isospin doublet of (non-relativistic) nucleon helds, 
isospin triplet of pion helds, and electromagnetic vector held are denoted by A, tt, and 
A^, respectively, and a and r are spin and isospin Pauli matrices. The arrow over the 
gradient specihes whether it acts on the left or right nucleon held. The isospin operators 
ejv and are dehned as 

Cat = (1-F r7/2 , kn = {ks + kvt^)/2 , = cn + , (3.20) 

and Ks and Ky are the isoscalar and isovector combinations of the anomalous magnetic 
moments of the proton and neutron. The power counting of the resulting vertices follows 
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by noting that each gradient brings in a factor of Q, so, for example, the two terms in 
//jtTV are each of order Q, while (ignoring the counting Q assumed for the external field 
A^) the hrst term in is of order Q, and the remaining ones in the second line of 

equation (13.1911 are of order 

In addition to the chiral Hamiltonians above, up to and including order (Q/A^^)^ 
there are fourteen contact interaction terms allowed by the symmetries of the strong 
interactions, each multiplied by a LEG. Two of these contact terms (proportional to the 
EEC’s Cs and Ct in standard notation) are of a non-derivative type, and therefore are 
of order (Q/A^)°, while the remaining twelve (proportional to the EEC’s C') of order 
(Q/A^)^ involve two gradients acting on the nucleon fields—they are listed in Ref. [143j . 
The contact potential at order [Q/Ayff' derived from them in the two-nucleon center- 
of-mass system in fact depends on Cs and Ct, and seven linear combinations of the 
C', customarily denoted as Ci,..., C 7 . The remaining Eve linear combinations of C[ 
have been shown to be related to C 5 and Ct by requiring that the Poincare covariance 
of the theory be satisfied to order (Q/A^)^ |143] . So the (Q/A^)^ potential involves 
nine independent EEC’s. (As a side remark, we note that the contact potential at 
order (Q/A^)'^ requires an additional hfteen independent EEC’s.) These EEC’s are 
determined by htting two-nucleon elastic scattering data. Minimal substitution in the 
gradient terms leads to a (two-nucleon) contact current |136l 1138] . 

Eastly, non-minimal couplings through the electromagnetic field tensor are also 
allowed. It can be shown |136] that the only two independent operator structures at 
order (Q/A^)^ are 


Tfirira 

ricT 


= e 


dx 


C[^N^(tN N^N + C\ 
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<TT 


•AT N^N 


- N^ctN 


■ V X A , 


(3.21) 


where the isoscalar C[^ and isovector C(g EEC’s (as well as the d'fs multiplying the higher 
order terms in the yvrA Hamiltonian) can be determined by fitting photo-nuclear data 
in the few-nucleon systems. We will return to this issue below. 


3.2.2. From amplitudes to currents We present the expansion of the amplitudes T and 
T-j, for the processes NN —> NN and NNj —> NN based on TOPT. Terms in this 
expansion are conveniently represented by diagrams. Before discussing them in some 
detail, it is worthwhile to make some preliminary considerations. We distinguish between 
reducible diagrams (diagrams which involve a pure nucleonic intermediate state) and 
irreducible diagrams (diagrams which include both pionic and nucleonic intermediate 
states). The former are enhanced with respect to the corresponding irreducible 
contributions by a factor of Q for each pure nucleonic intermediate state. In the static 
limit—that is, in the limit in which m —?■ 00 or, equivalently, neglecting nucleon kinetic 
energies—reducible contributions are infrared-divergent. The prescription proposed by 
Weinberg |112[ 1113] 1114] to treat the latter is to define the nuclear potential (and 
currents) as given by the irreducible contributions only. Reducible contributions. 
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instead, are generated by solving the Lippmann-Schwinger (or Schrodinger) equation 
iteratively with the nuclear potential (and currents) arising from the irreducible 
amplitudes. 

The formalism developed by some of the present authors is based on this 
prescription. However, the omission of reducible contributions from the dehnition of 
nuclear operators needs to be dealt with care when the irreducible amplitudes are 
evaluated in the static limit approximation, which is usually the case. The iterative 
process will then generate only part of the reducible amplitude. The reducible part of 
the amplitude beyond the static limit approximation needs to be incorporated order 
by order—along with the irreducible amplitude—in the dehnition of nuclear operators. 
This scheme in combination with TOPT, which is best suited to separate the reducible 
content from the irreducible one, has been implemented in Refs. |136( I137( 1138] and 
is described below. The method leads to nuclear operators which are not uniquely 
dehned due to the non-uniqueness of the transition amplitude off-the-energy shell. The 
operators, while non unique, are unitarily equivalent, and therefore the description of 
physical systems is not affected by this ambiguity. We defer to Ref. |137] for a discussion 
of these technical aspects. 

Another approach, implemented to face the difficulties posed by the reducible 
amplitudes, has been introduced by Epelbaum and collaborators [118] . The method, 
referred to as the unitary transformation method, is based on TOPT and exploits 
the Okubo (unitary) transformation [144] to decouple the Fock space of pions and 
nucleons into two subspaces, one containing only pure nucleonic states and the other 
involving states which retain at least one pion. In this decoupled space, the amplitude 
does not involve enhanced contributions associated with the reducible diagrams. The 
subspaces are not-uniquely dehned, since it is always possible to perform additional 
unitary transformations onto them, with a consequent change in the formal dehnition 
of the resulting nuclear operators. This, of course, does not ahect the physical results. 

The two TOPT-based methods outlined above lead to formally equivalent operator 
structures for the nuclear potential and electromagnetic currents up to loop-corrections 
included |138] . which makes it plausible to conjecture that the two methods are closely 
related. However, this topic has not been investigated further. 

In what follows, we outline briehy the methods developed in Refs. |1361113711138] 
and sketch how nuclear operators are derived from transition amplitudes. We start from 
the conventional perturbative expansion of the NN scattering amplitude T, which reads 

oo X \ n—1 

(/|rK>H/|//.g( ^^-„^^„ g.) U). (3.22) 

where \ i) and |/) represent the initial and hnal NN states of energy Ei = Ef, EIq is the 
Hamiltonian describing free pions and nucleons, and Hi is the Hamiltonian describing 
interactions among these particles (see section [3. 2. R . The evaluation of this amplitude 
is carried out in practice by inserting complete sets of Hq eigenstates between successive 
terms of Hi. Power counting is then used to organize the expansion. 
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In the perturbative series, equation fl3.22p . a generic (reducible or irreducible) 
contribution is characterized by a certain number, say M, of vertices, each scaling 
as X (*=1, • • •, M), where a* is the power counting implied by the relevant 

interaction Hamiltonian and /9j is the number of pions in and/or out of the vertex, 
a corresponding M — 1 number of energy denominators, and possibly L loops. Out 
of these M — 1 energy denominators, of them will involve only nucleon kinetic 
energies, which scale as Q^, and the remaining M — — 1 will involve, in addition, 

pion energies, which are of order Q. Loops, on the other hand, contribute a factor 
each, since they imply integrations over intermediate three momenta. Hence the power 
counting associated with such a contribution is 

M 


He 

^1 = 1 

Clearly, each of the M 
1 
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— (A/—Afx — 1) —2Afx 


X Q 


3L 


(3.23) 


Mk — 1 energy denominators can be further expanded as 
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(3.24) 


where Ej denotes the kinetic energy of the intermediate two-nucleon state, the pion 
energy (or energies, as the case may be), and the ratio {Et — Ei)/oJt, is of order Q. The 
terms proportional to powers of {Ei — Ei)/uJt, lead to non-static corrections. 

The Q-scaling of the interaction vertices and the considerations above show that T 
admits the following expansion: 

^ _ ^Gmin') _j_ 25) 

where ~ Q"', and chiral symmetry ensures that Umin is hnite. In the case of the 
two-nucleon potential = 0. A two-nucleon potential v can then be derived which, 

when iterated in the Lippmann-Schwinger (LS) equation, 

V V Gq V V Gq V Gq V . , (3.26) 


leads to the on-the-energy-shell {Ei = Ef) T-matrix in equation fl3.25p . order by order 
in the power counting. In practice, this requirement can only be satished up to a given 
order n*, and the resulting potential, when inserted into the LS equation, will generate 
contributions of order n > n*, which do not match In equation (I3.26p . Go denotes 
the free two-nucleon propagator. Go = l/(T^i — Ej + irj), and we assume that 

w = + ^ ('3_27) 

where the yet to be determined is of order Q^. We also note that, generally, a 
term like Gou*-”^ is of order since Go is of order and the implicit loop 

integration brings in a factor Q^. Having established the above power counting, we 
obtain 


, (3.28) 

n(i) = - [n(o)Gon(o)] , (3.29) 

n(2)=T(2)- [^;(o)Gon(°)Gon(°)] 

_ [y(i)Gon(°)+n(°)GonW] 


(3.30) 
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The leading-order (LO) term, consists of (static) OPE and two (non-derivative) 
contact interactions, while the next-to-leading (NLO) term, is easily seen to 
vanish inz!, since the leading non-static corrections of order Q in to the (static) 
OPE amplitude add up to zero on the energy shell, while the remaining diagrams in T^) 
represent iterations of whose contributions are exactly canceled by Go 
(complete or partial cancellations of this type persist at higher n > 2 orders). The 
next-to-next-to-leading (N2LO) term, which follows from equation fl3.30p . contains 
two-pion-exchange (TPE) and contact (involving two gradients of the nucleon helds) 
interactions. 

The inclusion (in Erst order) of electromagnetic interactions in the perturbative 
expansion of equation fl3.22l) is in principle straightforward. The transition operator 
can be expanded as [137] 


= T^-3) + t^- 2) + 7d-i) + ... ^ (3.31) 

where is of order e (e is the electric charge). The nuclear charge, p, and current, j, 
operators follow from = ^4° p — A-j, where = (A^, A) is the electromagnetic vector 
held, and it is assumed that has a similar expansion as T^. The requirement that, 
in the context of the LS equation, matches Tj order by order in the power counting 
implies relations for the = A^ which can be found in Ref. |137] . similar 

to those derived above for v^'^\ the strong-interaction potential. 

The lowest order contributing to the charge operator has n = —3, 

p(-3) = e + (1 2). (3.32) 


There is no n = —3 contribution to j, and the lowest order (n 
single-nucleon convection and spin-magnetization currents, 

;(-2) _ ^ /oTV I ^ 'Tpz 


2m 


o TV 1 , • Ai" 

2 Ki - V tcTi X q — 


= —2) consists of the 
+ (1 2) , (3.33) 


where q is the momentum carried by the external held, kj and Kj denote hereafter the 
combinations of initial and hnal nucleon momenta 


ki = Pi-Pi, Ki = (p'+ p,)/2 . (3.34) 

Note that and are the same operators as the one-body terms used in the 

conventional approach, see equations (13.41) and (13.bh . The counting e Q~^ (e Q~^) in the 
charge (current) operator follows from the product of a factor e (e Q^) associated 
with the 'jNN vertex, and a factor Q~^ due to the momentum-conserving 5-function 
implicit in a disconnected term of this type. 

The contributions to the electromagnetic current and charge operators up to one 
loop are illustrated diagrammatically in hgures[3]and|ll As already noted, the LO starts 
at n = —2 for the current and at n = —3 for the charge operator; NnLO corrections to 
both of them are labelled as Q” x LO. 

The currents from LO, NILO, and N2LO terms and from N3LO loop corrections 
depend only on the known parameters qa and /tt (NILO and N3LO), and the nucleon 
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LO: 


NILO: 


N2LO: 


N3LO: 


(a) 



(b) (c) 


Jf 

(d) 



(k) (I) (m) 


(n) (o) 


Figure 3. Diagrams illustrating one- and two-body currents entering at eQ~^ (LO), 
eQ~^ (NILO), eQ° (N2LO), and eQ^ (N3LO). Nucleons, pions, and photons are 
denoted by solid, dashed, and wavy lines, respectively. The square in panel (d) 
represents the (Q/m)^ relativistic correction to the LO one-body current; the solid 
circle in panel (j) is associated with the 'jitN current coupling of order e Q, involving 
the leg’s dg, dg, and dgi; the solid circle in panel (k) denotes two-body contact terms 
of minimal and non-minimal nature, the latter involving the LEC’s and C[q. Only 
one among all possible time orderings is shown for the NLO and N3LO currents, so 
that all box contributions also include crossed box contributions. 


magnetic moments (LO and N2LO). Unknown LEC’s enter the N3LO OPE contribution 
involving the ■jttN vertex of order e from second line of equation (I3.19p . as well 

as the contact currents implied by non-minimal couplings, equation (I3.2ip . discussed in 
the next subsection. On the other hand, in the charge operator there are no unknown 
LEC’s up to one loop level, and OPE contributions, illustrated in panels (c)-(e) of 
hgure m only appear at N3LO. The contributions in panels (d) and (e) involve non¬ 
static corrections |137] . while the contribution in panel (c) is associated with the yvriV 
coupling of order e Q originating from the hrst term in equation fl3.19p . It leads to a 
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Figure 4. Diagrams illustrating one- and two-body charge operators entering at e Q~^ 
(LO), eQ~^ (N2LO), e (N3LO), eQ^ (N4LO). The square in panel (b) represents 
the {Q/mY relativistic correction to the LO one-body charge operator, whereas the 
solid circle in panel (c) is associated with a 771 TV charge coupling of order e Q. As in 
figure m only one among the possible time orderings is shown for the N3LO and N4LO 
charge operators. 


two-body charge operator given by 
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[Ti ■ T2 -b T2z) 
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(3.35) 


In the present x^FT context, was derived first by Phillips in 2003 |145] . However, 
the operator of equation fl3.35p is the same as the vr-exchange contribution derived 
within the conventional approach described in section 13.1.21 (see Ref. [6] and references 
therein). This operator plays an important role in yielding predictions for the A=2-4 
charge form factors that are in very good agreement with the experimental data at low 
and moderate values of the momentum transfer (Q < 5 fm“^) [6l I138j . The calculations 
in Ref. |138j also showed that the OPE contributions from panels (d) and (e) of figure |4] 
are typically an order of magnitude smaller than those generated by that in panel (c). 

The loop integrals in the N3LO and N4LO diagrams of figures [3] and 0] are ultraviolet 
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divergent and are regularized in dimensional regularization |136l I137j . In the current 
the divergent parts of these loop integrals are reabsorbed by the LEC’s C- [136] . In the 
charge, however, they cancel out, in line with fact that there are no counter-terms at 
N4LO [137j . Finally, the resulting renormalized operators have power-law behavior for 
large momenta, and must be further regularized before they can be sandwiched between 
nuclear wave functions. This is accomplished by the inclusion of a momentum-space 
cutoff of the type Cf^{k) = exp(—/c"^/A^) with A in the range ~ (500-700) MeV/c. The 
expectation is that observables, like the few-nucleon form factors at low momentum 
transfer of interest here, are fairly insensitive to variations of A in this range. 

3.2.3. Determining the LEC’s We now turn our attention to the determination of the 
LEC’s. The Ci in the minimal contact current, corresponding to the A cutoffs of 500 
and 600 MeV/c in C^ik), are taken from hts to the two-nucleon scattering data |120j . 
The d/, entering the OPE N3LO current, could be htted to pion photo-production 
data on a single nucleon, or related to hadronic coupling constants by resonance 
saturation arguments. Both procedures have drawbacks. While the former achieves 
consistency with the single-nucleon sector, it nevertheless relies on single-nucleon data 
involving photon energies much higher than those relevant to the threshold processes 
under consideration and real (in contrast to virtual) pions. The second procedure is 
questionable because of poor knowledge of some of the hadronic couplings, such as 
gpNN- Alternative strategies have been investigated for determining the LEC’s d[ as 
well as C[^ and |138] . In this respect, it is convenient to dehne the adimensional 
LEC’s df’^ (in units of the cutoff A) related to the original set via 

C'a = df/A- , 4 = 

C;, = </A‘, d’,^d^/Ar ci'l=4/A^ (3.36) 

s 

where the superscript S' or E on the d^ ’ labels the isospin of the associated operator. 

The isoscalar df and df have been hxed by reproducing the experimental deuteron 
magnetic moment /i^ and isoscalar combination of the trinucleon magnetic moments. 
Their values are listed in table [H The LEC df multiplying the contact current is rather 
large, but not unreasonably large, while the LEC df is quite small |138] . 

Table 1. Adimensional values of the isoscalar LEC’s corresponding to cutoff 

parameters A in the range 500-600 MeV/c obtained for the N3LO/N2LO Hamiltonian. 


A 

df 

df X 10 

500 

4.072 

2.190 

600 

11.38 

3.231 


The isovector LEC’s df, df, and df have been determined in the following three 
different ways (denoted as I, II, and III in table [2]). In all cases I-III, we have assumed 
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Table 2. Adimensional values of the isovector LEC’s corresponding to cutoff 
parameters A in the range 500-600 MeV/c obtained for the N3LO/N2LO Hamiltonian. 
In sets II and III the values of dtf have been fixed assuming a A dominance model. 
See text for further explanation. 


A 

<(i) 

dpi) 

<(11) 

<(11) 

dHlII) 

4(III) 

500 

10.36 

17.42 

-13.30 

3.458 

-7.981 

3.458 

600 

41.84 

33.14 

-22.31 

4.980 

-11.69 

4.980 


d'f/d^ = 1/4 as suggested by A dominance in a resonance saturation picture of the 
N3LO OPE current of panel (j) in figure [3l In set I, d^ and d^ have been constrained 
to reproduce the experimental values of the np radiative capture cross section anp at 
thermal neutron energies and the isovector combination of the trinucleon magnetic 
moments. This, however, leads to unreasonably large values for both LEC’s, and is 
clearly unacceptable |138] . In sets II and III, the LEG d^ is hxed by assuming A 
dominance (see equation (4.2) of Ref. jl38j L while the LEG df multiplying the contact 
current is htted to reproduce either anp in set II or in set III. Both alternatives 
still lead to somewhat large values for this LEG, but the degree of unnaturalness is 
tolerable in this case. There are no three-body currents at N3LO |136j . and therefore it 
is reasonable to £x the strength of the two-nucleon contact operators by fitting a three- 
nucleon observable such as and ijX. Note that the values of the LEC’s in table [U 
and [21 as well as the results presented in section 01 have been obtained using the chiral 
NN potential derived up to N3LO by Entem and Machleidt in Ref. m, augmented, 
in the A = 3,4 cases, by the three-nucleon interaction derived up to N2LO, in the local 
version of Ref. UZB. 

Finally, we conclude by noting that hadronic electromagnetic form factors need to 
be included in the nuclear yEFT charge and current operators. The latter could be 
consistently calculated in chiral perturbation theory (yPT) |146j . but the convergence 
of these calculations in powers of the momentum transfer appears to be rather poor. 
For this reason, in the results reported below for the few-nucleon form factors, they are 
taken from hts to available electron scattering data [138j . 

3.3. The covariant spectator theory (GST) 

The covariant spectator theory (GST) [31 1147[ [ITS] , when applied to few-nucleon form 
factors, starts from the premise that they can be calculated from a covariant held theory 
in which nucleons and the lightest mesons are treated as the effective degrees of freedom. 
A covariant equation is constructed that provides an approximate solution to the held 
theory, and from this the currents are determined and the form factors calculated. 

































Electromagnetic Structure of Few-Nucleon Ground States 


31 


(A) (B) 

Figure 5. (Color online) Panel A shows the two-body CST-BS equation where the 
amplitude for both nucleons off-shell in the final state can be calculated from the CST 
amplitude. Panel B gives a diagrammatic representation of the two-body CST (or 
Gross) equation for the NN scattering amplitude. In both panels, crosses indicate 
on-shell particles and the (green) solid box is the kernel. 




Figure 6. (Color online) Di¬ 
agrammatic representation of the 
antisymmetrized two-body one bo¬ 
son exchange kernel. 


3.3.1. CST two-body equations Figure E] shows a diagrammatic representations of two 
forms (discussed below) of the CST equation for two-nucleon scattering. All diagrams 
are Feynman diagrams, so iterating the CST equation gives the sum where 

M" is the Feynman diagram with n kernels connected by a two-nucleon propagator 
with one nucleon on mass shell. The CST equation differs from the more familiar 
Bethe-Salpeter (BS) equation |149] in that one of the intermediate nucleons (labeled by 
an X in the diagrams) is always on-shell, while in the BS equation both are off-shell. 
It is important to realize that both the BS and the CST equations are exact, provided 
the kernel is exact, so that the approximation to the dynamics lies in the choice of 
kernel, not in the choice of equation. In this way the CST is similar to the conventional 
approach described above; the kernel plays the role of a generalized potential, with a 
theoretical structure based on physical insight and parameters adjusted to £t the data. 
However, a signihcant difference here is that the structure of the kernel corresponds 
to a selection of Feynman diagrams from which the electromagnetic current can be 
determined consistently (as in yEFT). 

In the following, we use the notation 

dS'k m r _ f dr'ki m f d^k2 m 

A,., ^ J (^^ J (^^ ’ ^ 

where m is the nucleon mass and Ek = \/ im? -|- the on-shell energy. The specific form 
of the CST equation for the NN scattering amplitude M, with particle 1 on-shell in the 
initial state and both particles off-shell in the hnal state, as illustrated in panel (A) of 
£gure[5l is 

Mi2(pt,K;P)=Fi2(pt,p;;P)- [ V,2{pl,k,;P)S2ik2)M,2{k,,p[;P), (3.38) 

Jki 

where the momentum of the internal particle 1 is on-shell (so that ki = {F'fc,ki}), 
pI is off-shell (so that pi = {po^Pi})) and V is the antisymmetrized kernel. In the 
















Electromagnetic Structure of Few-Nucleon Ground States 


32 


rest frame P = Pi + P 2 = {W,0}, so the momentum of the off-shell particle 2 is 
k 2 = P — ki = {W — Ek, —ki}, and its propagator is S' 2 (/c 2 ) = — ^2 ~ [In many 

references, the relative momenta p = |(pi — P 2 ) are used in place of pi (and similarly 
for k), and sometimes G is used in place of S; use care in reading the literature.] Since 
pI is off-shell in equation (13.381) . this is referred to as the CST-BS equation; it gives an 
expression for the fully off-shell scattering amplitude once the CSX amplitude is known. 
To hnd the CSX amplitude, one merely sets pi pi, giving a closed equation for 

M [T^rr^ 

Mi2(pi,K;P)=Fi2(pi,K;P)- [ Vi2{pi,kp,P)S2{k2)Mu{kup[-,P). (3.39) 

Jki 

This is the equation illustrated in panel (B) of hgure [5l Since particle 1 is on shell 
throughout, M 12 may be dehned to be the matrix element of the Feynman scattering 
amplitude M. between positive-energy Dirac spinors of particle 1 


Ml2iPl,p'i;P) = Mx\',aa'{Pl,P'l]P) = Ug{p>i,X)Mfig',aa'{Pl,p'l]P)UfifTp\,X') . 


(3.40) 


The index 2 still refers collectively to the Dirac index of particle 2, {aa'}, but index 
1 may refer either to the helicities of particle 1 {AA'} or the Dirac indices {PP'} [the 
replacement of a Dirac index (a*, /S*, ...) by a helicity index (A*) always indicates 
a corresponding contraction with a positive-energy helicity spinor]. For applications, 
it is sufficient to also place particle 2 on mass shell in the initial state. The two- 
body interaction kernel V that enters equations (I3.38p or (13.391) . sometimes called 
the “potential” because of its close connection with the nonrelativistic potential, is 
constructed by explicitly antisymmetrizing the kernel V, as illustrated in hgure El In 
its Dirac form it is 

Vgf}'^aa'{Pl,ki;P) = ^ [Vgg',aafPi,ki;P) + r]iVag',ga'{P2,kuP)] , (3.41) 

where the factor pj = C(—)^’''^ (with J=0 or 1 the isospin of the NN state) and C = 1 
for bosons and —1 for fermions. 

For the calculation of the deuteron form factors, a relativistic description of the 
deuteron bound state is needed. The deuteron appears as a pole in the scattering 
amplitude M at P^ = M|. At the pole, the scattering amplitude can be written 


etHPi,p)eyAp[,p) 


Mxy,aa'{pl,p'i;P) = - 




Mj 


p2 


R\\' ,aa' {pi,Pi,P), 


(3.42) 


where G^‘^{pi,P) = G^{pi, is the vertex function for an incoming deuteron with 

four-momentum P and helicity Xd (including the charge conjugation matrix; for details 
see Ref. |152] ). and i? is a non-singular remainder. At the pole, these vertex functions 
satisfy a homogeneous integral equation of the type (I3.38p or (13.391) depending on 
whether or not both particles are off-shell. For the discussion of the form factors, it 
is useful to define a relativistic wave function, 4/, equal to 


'>lUh,P) = S„Ak2)5^J,A'‘uP)' 


( 3 . 43 ) 
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Figure 7. Three-body diagrams (with internal spectators on-shell) that contribute to 
the infinite class of successive two-body scatterings (represented by the ovals). These 
are not time-ordered diagrams, but as in all other diagrams in this section, the order 
of the interactions develops from right to left. In this example, all of the diagrams 
but the first contribute to the subamplitude where particle 1 (on the top) is the 
last spectator and particle 3 (on the bottom) is the first spectator. The first diagram 
contributes 


3.3.2. GST three-body equation In the absence of irreducible three-body interactions 
(which is the case for CST one-boson-exchange models), any three-body scattering 
amplitude T can be viewed as a sum of successive two-body scattering processes, with 
one particle uniquely identihed as the spectator (see £gure[7j). It can then be written 
T = Ylij=i ^*'^ 5 where each component is the total amplitude of all processes in which 
particles i and j are the spectators in the hnal and initial states, respectively. Summing 
these contributions leads to a coupled set of integral equations with a Faddeev-like 
structure, but because these equations sum Feynman diagrams and not the time-ordered 
diagrams used in Hamiltonian theories, their physical content is quite different |153] . 

As in the two-body case, a three-body bound state with mass Mt appears as a pole 
in the scattering amplitude T aX = M|, where the conserved total four-momentum 
is P = ki + k 2 + ^ 3 , with ki the individual particle momenta. As in equation fl3.42p . 
near the pole the subamplitudes can be written 

|P)(F^’| 


jnj _ 


Mi - P2 




(3.44) 


where F* is a vertex function with particle i as spectator, and is a non-singular 
remainder. These vertex functions F* satisfy a homogeneous integral equation of the 
Faddeev type, whose kernel contains the total two-body scattering amplitude [5] . 

In the CST, the spectator is always on mass shell. Because the structure of the 
Faddeev equation dictates that after each two-body interaction a different particle 
becomes spectator, a second particle j must be placed on-shell in order to obtain a 
closed set of equations for the CST vertex function. We will use the notation F* to 
denote the subamplitude where particle i is the spectator and j ^ i is the second on- 
shell particle (which must be one of the particles in the “last” interacting pair). When 
the three particles are identical, the two-body amplitudes in the kernel are symmetrised, 
and the various vertex functions F* are related by symmetry relations. As a consequence. 
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Figure 8. (Color online) Diagram¬ 

matic representation of the Faddeev equa¬ 
tion (I3.451) . The interior dot labels the spec¬ 
tator to the “last” interaction. 


it is sufficient to determine just one of them, by convention (the indices 1 and 2 on 
r will be suppressed from here on). This vertex function describes the structure of a 
three-body bound state in the CST, and it also determines the three-body form factors. 

For three identical nucleons with mass m, the bound-state Faddeev vertex function 
to be calculated is Tx^x^aih, h, h) = Uaiih, >^i)ua 2 {k 2 , >^ 2 )'^aia 2 a{ki, k 2 , h), where 
nucleons 1 and 2 are on mass shell {kf = = m^), with helicities Ai and A 2 , and 

particle 3 is off mass shell with the associated Dirac index a. In the following the 
alternative notation F(/ci,/C 2 ; D) is used for the vertex function; it is more convenient 
because k^ = P — ki — k 2 is a. dependent variable. 

The vertex function satishes the Faddeev equation (shown diagrammatically in 
hgure [8]) 

rAiA2a(^l) ^ 2 ;-P) = — ^ [ AFA2A^,aa'(^2, ^ 2 ;-P 23 ) 2 CP ’12 fl^AiA^a'(^l) ^ 2 ;-P) ) (3.45) 

A' >^"2 

where P 23 = P —A;i is the total four-momentum of the pair, M is the two-body scattering 
amplitude satisfying equation fl3.39p . ^ = -|-1 for bosons and C = —1 for fermions, and 
P 12 is the permutation operator that interchanges particles 1 and 2 (for details, see 
Ref. 0 )- As in the two-body case, it is useful to define a three-body relativistic wave 
function T as 

^AiA2q(^15 ^2l P) = Saa'{.kfi)Vx^x^a'iki, k2] P). (3.46) 

3.3.3. Properties of OBE models used with the CST equations It has been found that 
a one-boson-exchange (OBE) kernel, illustrated in hgure [6l can provide a precision £t 
to the data. The bosons required include the familiar six bosons vr, 77 , cxo, ui, p, a; plus, 
in some cases, the heavier axial vector isoscalar hi and isovector oi mesons. There 


Table 3. Mathematical forms of the bNN vertex functions, with 0(p) = {m — fi)l2m. 
The vector propagator is = g^iz—q^iq^/rnf with the boson momentum q = pi—ki = 
k2 -P2. 


r{h) 

eb 

Ai $§ A 2 

A(p, k) or A^(p, k) 

0+(s) 

— 

A 1 A 2 

Qs - Vs [0(p) -b 0(fc)] 

0-(p) 

+ 

A 1 A 2 

9pl^ -9p{^ - K) [0(P)7^ + 3^Q{k)] 

l-(u) 

+ 


9v [Y + - k)u] +9vVv [0(p)7^ + 3^Q{k)] 

l+(a) 

+ 

KKdpu 

9a'l^Y 
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Table 4. Values of the 13 parameters for the 6 bosons of Model IIB and the 16 
parameters for the 6 bosons of Model W16. All masses and energies are in MeV; other 
couplings are dimensionless; Gb = gllAir. Parameters in bold were varied during the 
fit; those labeled with an * were constrained to equal the one above. The triton binding 
energy is Et (with its experimental value in parentheses). For model W16, Vs = 

No Et was calculated for model IIB. 




Model IIB (1992) 



Model W16 (1997) 


b 

I 

Gb 

mb 



Gb 

mb 

Xb/r^ Kv 


TT®, TT^ 

1 

13.37758 

138.0 

— 

— 

13.34 

138.0 

0.0 — 

2075 

V 

0 

5.30321 

548.8 

— 

— 

2.969 

548.8 

0.0 — 

1206 

ao 

0 

4.86870 

522.0 

— 

— 

4.99887 

506 

-1.2 

1206* 

Cl 

1 

0.24372 

482.0 

— 

— 

0.62818 

512 

4.16 

1206* 

UJ 

0 

8.86086 

782.8 

1.0 

0.22069 

14.879 

782.8 

0.0 0.195 

1206* 

P 

1 

0.60318 

760.0 

0.82989 

5.66983 

0.899 

760.0 

1.556 6.267 

1206* 


Ajv = 1675; Aj = 2185 

Ajv = 

1822; 

Et = -8.49 (-8.48) 


Table 5. Values of the 27 parameters for WJC-1 with 7 bosons and 2 axial vector 
contact interactions, and the 15 parameters for WJC-2 with no axial vectors, obtained 
by fitting np data. For additional explanation, see table SI 




Model WJC-1 (2008) 



Model WJC-2 (2008) 


b 

I 

Gb 

mb 

Xb/vb 



Gb 

mb 


Ah 

7T° 

1 

14.608 

134.9766 

0.153 

— 

4400 

14.038 

134.9766 

0.0 — 

3661 

TT^ 

1 

13.703 

139.5702 

-0.312 

— 

4400* 

14.038* 

139.5702 

0.0 — 

3661* 

V 

0 

10.684 

604 

0.622 

— 

4400* 

4.386 

547.51 

0.0 — 

3661* 

Co 

0 

2.307 

429 

-15.169 

— 

1435 

4.486 

478 

-2.594 

3661* 

Cl 

1 

0.539 

515 

4.763 

— 

1435* 

0.477 

454 

9.875 — 

3661* 

UJ 

0 

3.456 

657 

0.843 

0.048 

1376 

8.711 

782.65 

0.0 0.0 

1591 

p 

1 

0.327 

787 

-1.263 

6.536 

1376* 

0.626 

775.50 

-2.787 5.099 

1591* 

hi 

0 

0.0026 

— 

— 

— 

1376* 

— 




ai 

1 

-0.436 

— 

— 

— 

1376* 

— 





A 

jv = 1656; 

II 

1 

GO 

00 

(-8.48) 


Ajv 

= 1739; 

Et = -8.50 (-8.48) 



is a long and continuing debate as to whether the scalar mesons aj (where I is the 
isospin) are “real” mesons, but there is certainly a strong two-pion interaction in the 
isoscalar channel and it is well known that the nuclear force requires a signihcant two- 
pion exchange mechanism beyond that obtained from the iteration of two noninteracting 
pions. In the CST these mesons are treated as important phenomenological degrees of 
freedom with an effective mass chosen to £t the data (which comes out in the vicinity 
of 400-500 MeV). The simplicity of OBE models makes them very attractive, and they 
have been perused for over 40 years, but only with the recent fits using the CST equation 
has it been found that such models can give precision fits to the data. This success may 
be due in part to the cancellation theorem |154] . discussed briefly in section l3^ below. 
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The OBE kernels are the sum of one-boson-exchange terms of the form 

ki) (g) Al{p2, /C2) 


Vf,ipi,kv,P) = e,6- 


m‘, 


+ \q^ 


-f{^b,q) , 


(3.47) 


with b = {s,p,v,a} denoting the boson type (scalar, pseudoscalar, vector, and axial 
vector), q the four-momentum transfer, the boson mass, ej, a phase factor, 5 = 1 
for isoscalar bosons and 5 = Ti ■ T 2 = —1 — 2(—)^ for isovector bosons. The use of 
the absolute value |g^| amounts to a covariant redehnition of the propagators and form 
factors in the region > 0, done to remove all singularities [151] . The boson form 
factors /(Aft, q) are normalized to unity at = 0 and approach zero at large g^ with a 
mass scale of Aft, treated as one of the OBE parameters. The axial vector bosons are 
treated as contact interactions, with the structure as in equation 03.471) . but with the 
propagator replaced by a constant, where the nucleon mass m sets 

a convenient scale not related to a boson mass (the effective boson mass in a contact 
interaction is inhnite). The explicit forms of the numerator functions Aj 0 A 2 can be 
inferred from table El Note that Ap = 0 corresponds to pure pseudovector coupling, and 
that the dehnitions of the off-shell coupling parameters A or z/ differ for each boson. 

The bNN vertex functions, as written in equation fl3.47p . also include strong nucleon 
form factors 

A^(p, k) = h{p)h{k)A^{p, k) , (3.48) 

where p {k) is the four-momentum of the incoming (outgoing) nucleon, h{p) is a strong 
form factor associated with an off-shell nucleon with four-momentum p, and A^ is a 
reduced bNN vertex containing no form factors. The strong form factors h provide both 
needed convergence and a phenomenological description of the short range structure of 
the off-shell nucleons. The recent models |151] use the form 

(A^ -- 


h{p) = 


JL 

(A^ — m^)^ -|- (m^ — p^ 

with h{p) a function of p^ only, subject to the constraint that h{p) = 1 when = m^. 

The three types of OBE models used in applications of the GST equations are 
summarized in tables Eland El The oldest models (1992) |150] did not have any off-shell 
terms in the scalar meson exchanges (all z/ = 0) and Model IIB (one of four models 
developed at that time) ht the np data with a x^/Adata = 3.40, not a precision ht 
by todays standards (a revision of this model with slightly altered parameters gave a 
better ht with x^/Afdata = 2.53 |155] ). Shortly afterward, in connection with the study of 
the three-nucleon bound state |156] . it was realized that including the z^-dependent off- 
shell terms in the scalar exchanges would greatly improve both the ht to the two-body 
data and the three-body binding energy. Model W16 gave both the best ht to the np 
scattering data (with x^/Ndatm = 2.25) and the correct three-body binding energy |157j . 

The parameters for both of these older models were determined by hts to the 
Nijmegen phase shifts [158] , and it was found that htting these phase shifts did not give 
the best ht to the data itself. A major ehort was made to ht the np data base (below 


m 




(3.49) 
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Figure 9. (Color online) Results for calculations of x^/-^data (curves with a minimum 
and left scale) to a 2007 np data base, and triton binding energy Et (linear curves and 
right scale) for WJC-1 family (left panel) and for WJC-2 family of models (right panel) 
as function of i/a- The points with the lowest x^/iVdata are models WJC-1 and WJC-2, 
respectively. For further details, see Ref. |151j . 


Eiah = 350 MeV) directly, which produced the two precision hts (with x^/Ndata — 1) 
WJC-1 and WJC-2 nHi. Model WJC-1 was constructed with the goal to obtain the best 
possible ht, while the objective of WJC-2 was to use the smallest number of parameters 
without destroying the quality of the ht signihcantly. The variation of the three-body 
binding energy and the /Ndata for both of these models is quite striking and is shown 
in hgure [9l As with the family of models associated with W16, the best ht to the 
two-body data also gave the correct three-body binding energy. This result is a simple 
demonstration that the CST OBE models, which generate no irreducible three-body 
forces when applied consistently to the three-nucleon sector, seem to automatically 
include the physics contained in the explicit three-body forces required by nonrelativistic 
models 1I53|. 


3.3.4- Two-body current operator and the deuteron form factors The basic idea of the 
derivation of the current is to couple the photon to all propagators and momentum 
dependent couplings in each of the inhnite series of two-nucleon diagrams. This yields, 
according to a general argument developed by Feynman, a conserved current. If all 
contributions before and after the interaction with the photon are summed up using 
the two-body CST equation, hgure [5l it emerges that only the four diagrams shown in 
hgure [ini are needed to fully describe elastic electron scatting from the two-body bound 
state |15911160] . These can be written 

K.(k,Pp\jUP+.P-)h..KS(^-k') + VX,„.(kP^-,k’Pp]'ti,.X,P-) 


'kk> 


P+) (fc_, P_) k(k+, K) k.) S,,y{k. 


E (MeV) 
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(A) (B_) (C) 


Figure 10. (Color online) Diagrammatic representation of the current operator of 
the CST with particle 1 on-shell. 


n) P-) SAkAj^Af^+, k.) (k_, a;) (3.50) 


Jk- 

where here, and in the equations to follow, sums over repeated polarization indices 
are implied, and is the two-nucleon interaction current. The (B±) diagrams need 
the vertex function QagikA) for both particles off-shell (the CST-BS vertex function). 
This vertex function allows ^ m^, but when kf = it is related to the vertex 
function of equation (I3.42p by Qaxik, P) = Q^gik, PAgik, A). It can be calculated from 
the CST-BS equation fl3.38lh represented diagrammatically in panel (A) of figured As 
suggested by the hgure, the OBE kernels, obtained from Feynman diagrams, are already 
defined off-shell, and therefore once the CST vertex function is known, the CST-BS 
vertex function can be calculated by quadratures. These off-shell vertex functions are 
necessary to properly account for the interactions of the photon with particle 1, where 
energy-momentum conservation inside the loop does not allow the on-shell constraint 
to be maintained simultaneously before and after the interaction. Since the dnp vertex 
function is explicitly antisymmetric, it might be possible to replace these diagrams 
by (A)-type diagrams, but this would require an unconventional redefinition of the 
interaction current contribution (C) and has not been attempted. 

As shown by Riska and Gross (RG) [159] . the two-body current operator dehned in 
figure [lO] (and its generalization to inelastic scattering) will yield a conserved current, 
even in the presence of phenomenological form factors, provided the one- and two-body 
nucleon currents satisfy the appropriate Ward-Takahashi (WT) identities. In order for 
the RG prescription to work, the strong form factors at the meson-baryon vertices must 
be factorizable into a product of individual form factors associated with each hadron, 
and the models discussed above all have this property. 

With this structure, the strong nucleon form factors may be moved from the kernels 
to the nucleon propagators, leading to dressed (or damped) propagators of the form 

1 S-\p) _ 2m0(p) 

h2(p) h2(p) h2(p) 



where the h occurs squared because one comes from the initial and one from the final 
interactions that connect the propagator and 0(p) = {m — p — A)/(2m). A conserved 
two-nucleon current can then be constructed using a dressed single nucleon current of 
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the form [16111152j 

f{p,p) = h{p)h{p')jUP^p') ■ 

The reduced current satisfies the WT identity 
q^fniP^p') = eo\Sf\p') - Sf\p) 


(3.52) 


(3.53) 


where cq = e(l + Tz )/2 is the charge operator. 

The simplest solution to equation fl3.53l) . for the off-shell isoscalar part of the 
nucleon current needed for the description of the deuteron form factors, can be written 


as 


u r 


j'"ip',p) = /o(p',p)^ + £/o(p',p)0(p')^0(p) 


+ eofo{p',p) < FiY + F2 


2 m 


eo9o{p',p)F3Q{p')YQ{p), (3.54) 


where now cq = |e, g = p' —p,Fi = Fi[q^) are the isoscalar form factors of the nucleon, 
with F 3 a new form factor that contributes only when both nucleons are off-shell, and the 
transverse gamma matrix is Using the shorthand notation h = h{p) 

and h' = h(p'), the functions /o and go are 


k{p\p) = 


h' {mS — pk h {mS — p 
h p'"^ — p‘^ ^ h' — p’' 


9o{p\p) = 


AmS 


p/2 _ p2 


h 

hi 


h' 

X 


(3.55) 


The apparent singularity at g^ = 0 in the first line of equation (I3.54p is cancelled by 
the second line, because both Fi and F 3 are subject to the constraints Ti( 0 ) = 1 (for 
i = 1,3). Note that equation fl3.54p displays the very interesting property that all the 
physics is in the transverse term (second line), and that the longitudinal part (first line), 
uniquely fixed by the WT identity, does not contribute to a physical amplitude, because 
qU gives zero when contracted into another conserved current or into any of the real or 
virtual photon polarization vectors. 

The condition that the reduced interaction current U^(/cP+; k'PJ) (the current with 
the strong form factors h removed) must satisfy in order that the total current, J^{q), 
be conserved (the two-body WT identity) can be written as 


qu^g/ 3 ',aa'i^ P-) = eo Vg,aa> {k, F] P-) - Vgp'{k, F] P+) 

+ Vpfj'^aafk - g, F] P-) - Vggi^aafkF g; P 4 


(3.56) 


Recently the isoscalar part of this interaction current, needed for calculation of the 
deuteron form factors, has been uniquely determined |152] . 


3.3.5. Three-nucleon form factors in the CST To expose one of the central issues in 
the construction of three-body currents, we begin by looking at what appears to be the 
lowest order result in the BS formalism, and show that this expected result leads to 
overcounting. 

As discussed above, the full BS three-body vertex function is the sum of three 
subamplitudes T*, as shown in panel (A) of figure [HI Guided by nonrelativistic 
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Figure 11. (Color online) The figures enclosed in the rectangle show the BS vertex 
function as a sum of three subamplitudes for the three different choices of the last 
interacting pair. The rest of the figure shows that naive use of this vertex function in 
the form factor gives overlap terms (a)x(b)=(ab) and (b)x(a)=(ba) leading (after use 
of the equation) to two terms of type (d), where only one should be present. 


theory, we might expect the impulse approximation to the current to be related to 
the square of the wave function, as illustrated in panel (B). However, if this proposed 
current is expanded using the wave equations, it leads to two terms of type (d), while 
direct examination of the ladder sum (for example) shows that there should be only 
one such term. Unless an interaction term of type (d) is explicitly subtracted from 
the “impulse” approximation, it will be double counted. The same problem does 
not arise in nonrelativistic theory because there the diagrams represent a sequence of 
operators which, in general, do not commute. The iteration of (a) x (b) gives a different 
contribution from that of (b)x(a), and both must be present. (The treatment of this 
problem in the context of the BS theory is discussed in Refs. [16211163] .) 

It turns out that the spectator theory, like nonrelativistic theory, also does not 
suffer from double counting. Furthermore, the topology of the terms shown in hgure fTTl 
can be used to simplify the spectator formalism. A detailed demonstration of this is 
given in Refs. |164[ 1165] . 

In practice, the three-nucleon current is constructed by hrst coupling the photon 
to internal lines and vertices, which can be done very nicely using the “gauging of 
equations” method of Kvinikhidze and Blankleider [166] . In a second step, the three- 
body equations are employed to rearrange the result into a more usable form. For 
example, the hnal expression of the current obtained in Ref. [166] includes a contribution 
in which the spectator (particle 1) is off-shell. In Ref. [164] . the three-body bound-state 
and scattering equations are used to replace this amplitude by an equivalent one in 
which the spectator is on-shell, as reqnired in the GST. This replacement also leads to a 
nice demonstration of how the spectator equations avoid the double-counting problem 
in a natural way. The hnal result for the three-nucleon current is shown in hgure [T2j 
The diagrams (A)-(F) are referred to as the complete impulse approximation (CIA), 
while diagrams (G)-(J) are contribntions from the interaction currents. 
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Figure 12. (Color online) The electromagnetic three-nucleon current in CST for 
elastic electron scattering from the three-nucleon bound state. 


In algebraic form, the six diagrams that make up the CIA can then be written 

•^ciA = 3e / \^\iX2a'{ki,k2,kf)Oi2jl^,^{kf,k3)'^Xi\2a{ki,k2,k3) 

Jkik2 

+ ^Xi0'a{ki, kf, ks) Sg^fi{kf)j’^^{kf, fca) C>i2 u^{k2, Aa) 4 'aiA2«(^i, ^ 2 , ^ 3 ) 

+ ^XiX2a{ku k2, kf)u^{k2, A2) O12 kf) Sg>p{kf)Tx^i 3 a{ki, kf, , ( 3 . 57 ) 

where C>i 2 = 1 + 2 (Vi 2 j Ja'a^k'-ik) is the single-nucleon current discussed above, 
kf = ki±q, and in every term kf = k^ = ni^ and P = ki + k2 + k^i'S the four-momentum 
of the incoming three-nucleon bound state. The interaction current diagrams give 


= 3e f T 

Jkik2k^ 

where P' = P + q = ki + k'2 + k'.^, kf = k'f = k2 = rn'^, P^^ = k!^^ + k^\ and is the 
symmetrized two-body interaction current satisfying the two-body WT identity fl3.56p . 

Note that this calculation requires knowledge of the three-nucleon vertex function 
with the two interacting nucleons off-shell. This vertex function was defined in Ref. [164] 
and can be obtained using the Faddeev equation (I3.45p . generalized to the case when both 
of the hnal state interacting nucleons are off-shell. It requires the scattering amplitude 
for both of the final particles off-shell, which is obtained by quadratures from the off- 
shell kernel and the on-shell scattering amplitude using equation (I3.38p . illustrated in 
figure [5](A). The off-shell kernel is known (in principle), and is discussed in more detail 
in Refs. [T571IT^ fWf] . 

When the three-nucleon form factors were calculated with high-precision models 
WJC-1 and WJC-2 |168j . the three-nucleon vertex functions with two nucleons off-shell 
were not available, and were replaced by vertex functions with only one off-shell nucleon. 
This approximation was called “CIA-0”, because it can be interpreted as the zeroth- 
order Taylor expansion of the vertex function in the momentum of one nucleon around 


X'^a'{kl, k2, kf)Oi2Vy . ^2-^23)C>12^AiA2a(^1, ^2, ^ 3 ) ; (3.58) 
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its on-shell point. To test its quality, CIA-0 was compared with CIA for the case of the 
W16 model and found to be an excellent approximation [168j . 

When g —)■ 0, the nucleon propagators in the second and third terms of the CIA 
result of equation fl3.57p develop singularities that cancel, leading to terms involving the 
derivatives of the two-body kernel. Diagrams (B) and (C), as well as (E) and (F), of 
hgure [12] should therefore always be considered together to allow these cancellations to 
take place. The same situation occurs in the two-body case, where the (Bj-) diagrams 
of hgure HO] must be considered together. 

3.4- Appraisal of the different theoretical approaches 

The conventional and yEFT approaches use Hamiltonian dynamics, and treat relativity 
and electromagnetic gauge invariance perturbatively, expanding in powers of the ratio 
of the typical momentum of the process to the nucleon mass. Using yEFT, the 
potential is expanded in a power series. The advantage of these methods is that the 
underlying formalism is well-known and familiar, and the perturbative treatment of the 
potential permits an estimate of the theoretical error. The disadvantage is that the 
range of convergence of the yEFT series (the chiral symmetry breaking scale, A^ ~ 1 
GeV) is low for some of the calculations of few-body form factors presented here, 
which extend to momentum transfers beyond 1 GeV. The method in its purest form 
should produce results independent of the cutoff scale, because all divergences can be 
absorbed into parameters that encode the short range physics that cannot be calculated 
dynamically |169j . yet in the practical applications reviewed here, the requirement of 
cutoff independence limits the cutoff to a range of about 500 to 600 MeV, near the scale 
of the important physics described by two-pion exchange. Furthermore, the number of 
short-range parameters is not small; there are 24 undetermined parameters needed for 
the nuclear force at N3LO or [Q /and 5 more for the electromagnetic current to the 
orders discussed in this review. While these parameters can, in principle, be calculated 
from QCD, today they must be £t to experiment. However, this is the only approach 
that can currently be used to describe nuclei for mass numbers A > 4. 

The CST approach is almost completely complementary to the yEFT approach 
(which, for the purposes of this discussion, includes the material presented in section I3T] 
as well as section l3^ . It is based on covariant held theory and requires the development 
of a less familiar formalism (using covariant equations with relativistic normalization 
conditions) needed to treat the nonperturbative features of hadronic interactions at 
high energy. An advantage of the formalism is that it is fully covariant, exactly gauge 
invariant, and includes a four-current consistent, to all orders, with the dynamics. 
Disadvantages are that there is no systematic way to construct the kernel (he., potential) 
as there is in yEFT, making it difhcult to estimate the theoretical error, and the needed 
regularization is achieved by using hadronic form factors with cutoff masses that depend 
very sensitively on the £t. However, the OBE models used with this approach can be 
very efficient: one model has only 15 parameters and yet gives a £t to the np data of the 
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Figure 13. (Color online) Diagrammatic 
representation of the cancellation theorem 
discussed in the text. The heavy baryon is 
represented by the thick line. 


same quality as the best conventional models. Unfortunately, the necessary equations 
have not yet been extended beyond the A = ?> sector, and hts to the pp data have yet 
to be done. 

These two methods (yEFT and GST) are so completely different from start to hnish 
that it is difficult to compare them. Here we look at a few examples of how they might 
be describing the same physics in completely different language. 

One key to this understanding may be the cancellation theorem, illustrated 
diagrammatically to 4th order in hgure [13] |154] . In cases where a neutral meson is 
exchanged between a heavy baryon with mass M and another baryon of mass m, the 
sum of the Feynman ladder and crossed ladder diagrams is well approximated by the 
single box diagram with the heavy baryon on its mass shell. Briefly, this comes about 
because the evaluation of both of these diagrams reduces to the sum of contributions 
from nucleon and meson poles, and in the heavy mass limit the “large” meson pole 
contributions tend to cancel. This cancellation works to all orders, and in the limit 
when M —)■ cxo it is exact. Therefore, the success of the CST-OBE models may be 
due in part to the fact that GST sums these leading contributions to all orders (recall 
hgure |5]). If a charged meson is exchanged the cancellation is not exact, but it has 
been shown that for charged pion exchange the residue is well approximated by the 
exchange of scalar mesons with a distributed mass [1481 llTOj . giving another reason 
why scalar meson exchange is so important. In yEFT, the same physics is included 
order-by-order through (i) explicit calculations of crossed pion boxes and other time 
orderings (including contributions from so-called “stretched boxes”), and (ii) non-static 
corrections arising from the iteration of static approximations to lower order irreducible 
potentials (recall the discussion in section I3.2.2|) . The conclusion to be drawn from 
this discussion is that a simple comparison of the physics included in the two methods, 
including isolation of any differences that might exist, is difficult. 

A second key to the comparison is illustrated in hgure [T4j The GST off-shell OBE 
couplings will cancel neighboring propagators, making it possible to reinterpret the 
iteration of off-shell contributions to OBE exchanges as contact interactions, which can 
generate many non-OBE contributions in the two-body sector and effective three-body 
force contributions in the three-body sector. As suggested by the hgure, this can happen 
in an inhnite number of ways, showing how some of the triangle and three-body force 
diagrams needed in yEFT emerge from GST-OBE models, even though they do not 
appear to be included anywhere. This insight probably explains why GST-OBE can 
predict the three-body binding energy (recall hgure |9|) without adding GST three-body 
forces, while they are essential if one uses yEFT or the conventional approach. 

As has been understood since the hrst days of the GST (see, for example. Ref. [147] ). 
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Figure 14. (Color online) Boson-nucleon vertices with off-shell couplings (those that 
depend on the operator 0) can generate effective NN force diagrams beyond OBE, and 
effective thre-nucleon forces. In these examples the off-shell scalar a meson coupling 
generates an effective triangle contribution to the NN interaction [panel (A)] and a 
double triangle contribution [panel (B)]. In the three-nucleon sector, off-shell exchanges 
between two different nucleons can also generate effective three-nucleon forces [panels 
(C) and (D)]. Similar mechanisms arise from the pseudovector ttNN vertex. 


the antiparticle contributions from the CST off-shell propagators also make higher order 
contributions similar to those just described. In the rest frame of a two-body system, the 
propagator of the off-shell particle 2 may be decomposed into its positive and negative 
energy contributions 


G{k) = 


m- 




le 


m 

^ A 


M(k, A)M(k, A) 'y(—k, A)n(—k, A) 


2E.-W 


W 


(3.59) 


where W is the energy of the two-body system and u (v) are the positive (negative) 
energy Dirac spinors of the nucleon with three-momentum k. Using this decomposition, 
the deuteron bound state equation can be written in a coupled channel form, which in 
coordinate space becomes, approximately, 

{2Ey — Mci)i}^{r) = — {r)if^ {r) — V^~ {r)il)~ [r) 

—Mdi>~{r)= —V~^(r)il)^(r) — V (r)il)~(r). (3.60) 

Eliminating 'ip~ gives an effective potential for the positive energy sector 

U+-(r)U-+(r)- 


(2Uv - Md)^/;+(r) = - U++(r) + 




(3.61) 


The short-range piece proportional to V~^ V is positive dehnite in central channels, 
and when combined with a modest w-exchange term gives a credible explanation of 
the short-range repulsion needed for an explanation of the NN data. Perhaps more 
signihcantly, this simple argument demonstrates how negative-energy (antiparticle) 
channels generate higher order terms that could be compared to the yEFT expansion. 
As is obvious from the discussion in the last two paragraphs, the comparison of the 
dynamics included in CST-OBE with that included in yEFT is a challenging problem 
that has yet to be seriously studied. 
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All of these considerations also apply to the construction of the current (which, 
in the language of the yEFT or conventional approaches, is separated into charge and 
three-vector current operators). The current should be given by the theory, and in this 
respect yEFT provides clear guidance for its construction. Unfortunately, both yEFT 
and CST depart from their underlying field theories by using the measured structure of 
the nucleon to describe the jN N vertices, thereby giving up the goal of describing 
the few body systems entirely from first principles. Furthermore, the internal couplings 
of the photon to the pion (and, in CST, other mesons) are also not calculated from first 
principles. In the isoscalar channel needed for the deuteron form factors, coupling to 
mesons in flight does not contribute to the currents, so the isoscalar interaction currents 
are more tightly constrained. The freedom that might have existed because of the 
momentum dependence arising from the off-shell (i^-dependent terms) in the CST has 
recently been constrained |152] . However, there are always exotic interaction currents, 
such as the famous pvry and aojj contributions HU]. These effects are purely transverse, 
and hence cannot be constrained in CST by the WT identities that link the currents 
to the NN interaction. The two isoscalar contact interactions that arise in the yEFT 
expansions include these unknowns, but the unknown form factors associated with these 
terms are still a problem. Fortunately, the data show that these effects are small. 

We conclude this discussion with the observation that the interactions between 
hadrons can be described using either confined quark-gluon degrees of freedom or 
physically observable hadronic (i.e. colorless) degrees of freedom. Which choice is 
most efficient depends on the energy scale of the physics being studied. In this respect 
yEFT and CST are similar: CST uses nucleons and mesons with masses below 1 GeV, 
while yEFT uses nucleons and pions (only). The great advantage of yEFT is the 
organizational principle it provides: the precise definition of a perturbation series for 
the potential rooted in the approximate chiral symmetry of QCD. Unfortunately the 
scale at which this series diverges is probably less than 1 GeV. The organizational 
principle behind the CST is similar to that in vogue for many years: the exchange of 
the lightest mesons should account for the longest range (and hence largest) force, as 
suggested by the famous Yukawa potential proportional to 

-m(,r 

H(r) ~ Cb -. (3.62) 

r 

The longest range forces are described by the pion, and the intermediate and shorter 
range forces described by the exchange of heavier mesons (or, in the language of yEFT, 
multiple pion exchanges). This organizational principal is less precise and not as well 
defined as the one provided by yEFT, but works extremely well if CST equations are 
used with OBE parameters adjusted to fit the data. The most unifying conclusion to 
be drawn is that perhaps a comparison of how the physics is described by these two 
different approaches may teach us more than the study of either of them alone. This 
can be done by making a systematic nonrelativistic expansion of the CST equations, 
along the lines outlined in equations (I3.60p and (Id.blh : it is yet to be done. 
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4 . Results: theory vs. experiment 

In this section the theoretical predictions for the electromagnetic form factors of the 
Hydrogen and Helium isotopes corresponding to the conventional, yEFT, and CST 
approaches described in section |3] are compared to the experimental data obtained in 
the global analysis discussed in section El We first discuss the sensitivity of the various 
theoretical approaches to their physical content (section 14.ip . and then compare the 
predictions of the different approaches to data (section I4.2p . In section 14.21 we also 
present results for the charge and magnetic radii, and the magnetic dipole and electric 
quadrupole moments. 

4 . 1 . Model sensitivity 

Each of the theoretical approaches makes certain assumptions about the physics 
of electron scattering, and the final result is usually the sum of several different 
contributions. In this subsection we discuss the relative sizes of these contributions 
and the sensitivity of predictions to the values of unknown input parameters entering 
the various approaches. 

4-1.1. Conventional approach This approach uses the non-relativistic Argonne Uig 
(AVI8) two-nucleon potential jH2] for the A = 2 system. When applied to the 
A = 3 and 4 systems, it includes the (non-relativistic) Urbana IX (UIX) three-nucleon 
potential [92], whose strength parameters are adjusted to reproduce the ^H binding 
energy in exact Green’s function Monte Carlo calculations and the saturation density of 
nuclear matter in hypernetted-chain variational calculations. The resulting AV18/UIX 
Hamiltonian then leads to ^He and ^He binding energies in excellent agreement with 
the empirical values. Leading two- and three-body terms in the electromagnetic charge 
and current operators are constructed from the AV18 and UIX potentials, as outlined 
in subsection 13.11 and described in detail in Ref. [inO] . For the case of the deuteron only, 
we also present results obtained with an approximately relativistic treatment of nuclear 
dynamics, based on the relativistic Hamiltonian in equation (I3.14[l and by including one- 
and two-body terms in the electromagnetic operator as well as boost effects to order 
{v/cY in the initial and final states. 

Figure US] compares theoretical predictions for the deuteron form factors 
obtained with the AV18 potential using a one-body current only, and the full 
one plus two-body current operators. For this latter case results obtained with 
the approximately relativistic Hamiltonian are also presented. Two-body charge 
contributions, predominantly due to the pion-range operator of equation (13.hh . have 
opposite signs in Gc and Gq, substantially reducing Gc while moderately increasing 
Gq. 

The sensitivity of the three-body form factors to many-body terms in the 
electromagnetic current is displayed in figure [161 These many-body terms make 
important contributions for momenta larger than about 2 fm“^. 
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Figure 15. (Color online) The deuteron form factors calculated in the conventional 
approach are compared to data. The thin and thick (black) solid lines correspond, 
respectively, to results obtained by retaining in the charge and current operators only 
one-body (lA) or both one- and two-body contributions. The dash-double-dotted (red) 
line corresponds to the approximately relativistic calculation mentioned in the text, 
and includes contributions from one- and two-body electromagnetic operators. 


4-1.2. approach The xEFT calculations are based on the next-to-next-to-next- 

to-leading order (N3LO) two-nucleon potentials of Refs. |119[ 1120] corresponding to 
short-range cutoffs A = 500 and 600 MeV/c—denoted respectively as N3LO(500) and 
N3LO(600)—and retain, in the case of the A = 3 and 4 systems, the next-to-next-to- 
leading order (N2LO) three-nucleon potential in the local form of Ref. |121] . The low- 
energy constants (LEC’s) cd and ce (in standard notation) that characterize this three- 
nucleon potential have been constrained by reproducing, in essentially exact calculations 
based on hyperspherical-harmonics techniques |172] . the ^H/^He binding energies and 
the tritium Gamow-Teller matrix element for each of the A values considered—one of 
the LEC’s also enters the nuclear weak axial current and can therefore be determined 
in a weak transition. The resulting Hamiltonians are denoted as N3LO/N2LO(500) 
and N3LO/N2LO(600) below. Up to one loop, no unknown LEC’s enter in the 
electromagnetic charge operator, beyond the nucleon axial coupling constant qa, pion 
decay amplitude and proton and neutron magnetic moments, the latter associated 













Electromagnetic Structure of Few-Nucleon Ground States 


48 



Q [fm'‘] Q [fm ‘] 


Figure 16. ^He and charge and magnetic form factors as function of the momentum 
transfer Q (in fm“^), calculated within the conventional approach based on the 
AV18/UIX realistic potentials, are compared with experimental data. The thin and 
thick solid lines correspond to results obtained by retaining, respectively, one-body 
only or both one- and many-body contributions in the current and charge operators. 


with a next-to-next-to-leading order (N2LO) relativistic spin-orbit correction (see the 
N2LO panel in hgure Hj). In contrast, the electromagnetic current operator up to one 
loop is characterized by an additional hve EEC’s. In the yEFT results presented in the 
following subsections, the values listed in tables [Hand [2] (set III) are considered: the two 
EEC’s multiplying isoscalar operator structures are hxed by reproducing the deuteron 
and trinucleon isoscalar magnetic moments, while estimates for two of the three isovector 
EEC’s are obtained from A-isobar saturation arguments, with the remaining EEC hxed 
by reproducing the trinucleon isovector magnetic moment. 

Figure [T7] shows predictions for the deuteron form factors based on the N3EO(500) 
and N3EO(600) chiral potentials and including either EO electromagnetic operators only 
or corrections up to N3EO to these operators (see hgures [3] and |4]) . Eoop corrections at 
N4EO in the charge operator do not contribute to the observables under consideration, 
since they are isovector. The full results reproduce well the measured form factors 
for momentum transfers up to 2-3 fm“^. For Gq, the agreement between theory and 
experiment extends over a considerably wider region of momentum transfers, in fact 
well beyond the range that one would naively expect to be valid for the present yEFT 
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Figure 17. (Color online) The deuteron form factors calculated in the x^FT approach 
are compared to data. The thick and thin dashed (solid) lines are results for the 
N3LO(500) and N3LO(600) chiral potentials using LO (up to N3LO) corrections in 
the charge and current operators. 


expansion. The cntoff dependence is modest for Gq, but pronounced for Gq and Gm- 
Finally, we note that, as in the case of the conventional approach, the N3LO one- 
pion-exchange charge operator in panel (c) of hgure |4] gives the dominant contribution 
beyond LO. However, its effect on Gc and Gq is signihcantly smaller than obtained 
in the conventional approach. This is due to the fact that the functional forms of the 
cutoff used to regularize the operator in these two approaches are different and, more 
importantly, the value for A in the conventional AV18 calculation (A ~ 1200 MeV/c) is 
much larger than those adopted in the yEFT calculations (A either 500 or 600 MeV/c). 

Figure ITSl shows the sensitivity of the three-body form factors to the chiral cutoff 
and the order of the calculation. Here, since there are isovector contributions to the 
currents, corrections to the current up to N4LO are retained. 

4-1.3. GST approach Unfortunately, the most recent calculation for the deuteron form 
factors, reported in 1995 II73I, uses the older Model HB (with parameters given in 
table Hj). This model has no off-shell scalar meson couplings (the parameter Ug = 0, 
ensuring that the off-shell projection operators 0 given in table |3] are absent), yet these 
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Figure 18. (Color online) ^He and charge and magnetic form factors as function 
of the momentum transfer Q (in fm“^), calculated within the xEFT approach, 
are compared with experimental data. The thick and thin dashed (solid) lines 
correspond to results obtained using the N3LO/N2LO(500) and N3LO/N2LO(600) 
chiral potentials, with corrections in the charge and current operators up to LO (up 
to N4LO). 

terms are needed to give the high precision fits to the np scattering data |151] . Since 
the momentum dependence in these terms generates a new class of isoscalar interaction 
currents, a new generation of calculations that includes these currents is required. In 
2014, using principles of simplicity and picture independence, these isoscalar interaction 
currents were uniquely determined |152] . and excellent results for the deuteron magnetic 
moment |167] and quadrupole moment |174] were obtained (reported below). The new 
results for the form factors were not available at the time this review was prepared. 

Model IIB requires no interaction currents, but the results depend on the new, 
unknown, off-shell nucleon form factor, F 3 {Q), which contributes only when the 
incoming and outgoing nucleons are both off-shell. In has also been customary to add the 
contributions of a pvry interaction current, which is separately conserved and therefore 
not constrained by the requirements of current conservation. The size of these effects 
are shown in figure [191 where the the pvry form factor and coupling constant were taken 
from Ref. [175] (with the model 2 form factor). The hgure shows the extreme limits of 
FsiQ) = 0 or 1 [actually, F^i^Q) = 0 is impossible because of the constraint ^3(0) = 1, 
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Figure 19. (Color online) The deuteron form factors calculated in the CST approach 
are compared to data. The solid line (black) is the full calculation, the dashed line 
(red) has no exchange current (the CIA), the dash-dotted line (blue) is the extreme 
limit F3 = 1, and the dashed-dot-dot line (green) is the extreme limit = 0.) 

and is shown only to provide a lower limit, and while F^IQ) = 1 is possible, it violates 
our expectations that F 3 ^ 0 as Q —t 00 ]. Note that all of these effects are small except 
at high momentum transfers. 

The three body calculations are more recent |15711168] , and were done when the 
models WJC-1 and WJC-2 were available. Figure 1201 shows the isoscalar and isovector 
combinations of the A = ?> charge and magnetic form factors in the complete impulse 
approximation CIA-0 described in section 13.3.51 A moderate model dependence can 
be observed for momentum transfer above about 3 fm“^ (in the isovector magnetic 
form factor above 5 fm“^). The difference can be traced back to the different pion- 
nucleon coupling of the two models: the pseudoscalar admixture in WJC-1 automatically 
generates strong Z-diagram contributions suppressed by the pure pseudovector coupling 
used in WJC-2. 

To compare the CST calculations to the experimental data, at least the dominant 
pion-exchange currents would have to be taken into account, in particular the jirNN 
contact terms induced by pseudovector irNN coupling. They are part of diagrams 
(G) - (J) of hgure [121 which are not present in CIA. But Z-diagrams for pseudoscalar 
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Figure 20. (Color online) Isoscalar (left panels) and isovector (right panels) 
combinations of the charge (upper panels) and magnetic (lower panels) A = 3 form 
factors as function of the momentum transfer Q (in fm“^), calculated in the CST 
approach, are compared with experimental data. The results for the WJC-2 (red 
solid lines) and WJC-1 (blue dashed lines) were calculated in CIA-0. The results 
obtained within the conventional approach based on the AVI8/UIX potentials using 
the one-body operator and the same (Galster) nucleon form factors are also shown as 
dash-dotted lines. 


coupling are roughly equivalent to these pseudovector contact terms, so they are already 
partially included in the WJC-1 calculations in CIA. However, the corresponding mixing 
parameter has the opposite sign, so—compared to WJC-2—this contribution actually 
moves the form factors away from the data. In a complete calculation this will be 
compensated by a stronger contact term coming from WJC-l’s pseudovector nNN 
component. 

A direct comparison of a CIA calculation with data makes sense only for model 
WJC-2 in the case of the isoscalar magnetic form factor, where the 'ynMN contact term 
is suppressed. And indeed, as hgure EO] shows, the model describes the data very well. 
Figure l20] also compares the CST models in CIA with a calculation in the conventional 
approach performed with the AV18/UIX potential in impulse approximation with 
relativistic corrections and using the same (Galster) nucleon form factors. The close 
agreement between the results of AV18/UIX and WJC-2 can again be attributed to the 
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suppression of Z-diagrams in WJC-2. 

These results alone would not be sufficient to give preference to either WJC-1 or 
WJC-2. However, recent precision calculations of the deuteron quadrupole moment 
|174] seem to prefer WJC-2. For this reason, when comparing later on to the other 
approaches the H = 3 results with WJC-2 are shown. 

4-2. Comparison of the different approaches with data 

Electromagnetic form factors characterizing the coupling of the external held to 
individual hadrons enter the nuclear current and charge operators. Those for the proton 
and neutron in the dominant one-body terms of these operators are taken from hts 
to elastic electron scattering data off the proton and deuteron, specihcally the dipole 
parameterization (including the Galster factor for the neutron electric form factor) in 
the conventional approach, the Hohler parameterization |176] in the yEFT approach, 
and the GKex05 parameterization |177[ 1178] in the CST approach. We note that 
for momentum transfers up to ~ 6-7 fm~^, i.e., the Q-range over which most of the 
results are presented below, these various parameterizations do not differ signihcantly. 
Hadronic electromagnetic form factors also enter the nuclear many-body current and 
charge operators, and the specihc parameterizations adopted for these have been briefly 
discussed in section |3] and more extensively in the original references. 

4-2.1. nucleus Figure [21] compares the experimental data for the three deuteron 
form factors and for the A structure function (which has been measured out to 
~ 150 fm“^), to the full predictions obtained from the conventional, yEFT, and 
CST approaches. Note that the yEFT predictions for A only extend up to ~ 45 
fm“^, which is already well beyond the range of applicability of this approach. 

At low and moderate values of the momentum transfer (Q ^ 4 fm “^), the three 
theoretical approaches reproduce the data well. At larger values of Q, particularly in 
the diffraction region of the magnetic form factor and for the A structure function, 
the conventional and yEFT results are at variance with data, in the case of A 
by orders of magnitude at the highest Q’s. On the other hand, the CST results 
provide a remarkably good reproduction of the data over the whole Q-range covered 
by experiment—note that the A structure function drops by 9 orders of magnitude over 
this range! The quantitative success of the CST approach clearly demonstrates the need 
for a fully relativistic treatment of nuclear dynamics at high momentum transfers (Q ^ 5 
fm“^). It also suggests that, even in the extreme kinematical range covered by the A 
measurements, the description of nuclei in terms of protons and neutrons interacting 
via effective forces (and via effective currents with external electroweak helds) is far 
more robust than one would have naively expected. Indeed, there are no indications 
for the need to explicitly account for quarks and gluons, the degrees of freedom of 
the fundamental theory (QCD) governing their dynamics. It appears, instead, that 
the effects of the nucleon substructure can be subsumed in these effective forces and 
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Figure 21. (Color online) The deuteron experimental charge, quadrupole, and 
magnetic form factors compared to results obtained from (i) the conventional approach 
using the AV18 potential and including one- and two-body terms in the charge and 
current operators (thick solid line); (ii) an approximately relativistic version of the 
conventional calculation using equation (13.141) with the AV18 potential and including 
contributions from one- and two-body electromagnetic operators (dash-double-dotted 
line); (hi) the yEFT approach using rthe N3LO(500) (dashed line) and N3LO(600) 
(dash-dotted line) potentials and including corrections up to N3LO in the charge and 
current operators (with the region between these two lines lightly shaded); and (iv) the 
CST approach using model IIB (thin solid line). Also shown are conventional, ^EFT, 
and CST predictions for the A structure function measured out to ~ 150 fm“^ ~ 6 
GeV2. 


currents. 

The deuteron charge radius, and magnetic dipole and electric quadrupole moments 
are listed in table [6l The results in the conventional approach based on the AV18 
potential include one- and two-body terms in the electromagnetic charge and current 
operators, and under-predict the charge radius by 0.5%, the magnetic moment by 1.2%, 
and the quadrupole moment by 2.0%. We note that the isoscalar two-body current 
contributions from the momentum-dependent spin-orbit, L^, and quadratic spin-orbit 
components of the AV18 are individually small and tend to cancel out, since they have 
opposite signs. The pvry contribution is also found to be negligible. Indeed, the value for 
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fi{d) reported here, 0.847 n.m., coincides with that obtained in lA, i.e. with one-body 
currents only. The charge radius too is unaffected by isoscalar two-body contributions 
in the charge operator. However, these contributions increase the lA prediction for Q{d) 
by over 3%, but do not fully resolve the discrepancy between theory and experiment. 

Table 6. The deuteron charge radius rdd) (in fm), magnetic dipole moment fi{d) (in 
nuclear magnetons), and electric quadrupole moment, Q{d) (in fm^), obtained with the 
three different theoretical approaches, are compared to experimental values. The CST 
charge radius (marked with *) was calculated with the older model IIB, the magnetic 
and quadrupole moments with model WJC-2. The numbers in parentheses at the 
side of the WJC-2 calculations are an estimate of numerical errors. The experimental 
charge radius is taken from Ref. [12) . The experimental error for g{d) is not shown, 
since it is negligible. The numbers in parentheses at the side of the yEFT predictions 
for Q{d) and rdd) give the cutoff dependence of the results. The yEFT result for fi{d) 
is underlined, since g{d) is used to constrain one of the two EEC’s in the isoscalar 
current at N3LO. 


Observable 

Conv. 

yEFT 

CST 

Exp. 

rc{d) 

2.119 

2.126(4) 

2.085* 

2.130(10) 

/i(d) 

0.847 

0.8574 

0.864(2) 

0.8574 

Q{d) 

0.280 

0.2836(16) 

0.2836(3) 

0.2859(6) 


The yEFT approach leads to values for the deuteron static properties in very close 
agreement with experimental data. The value for the deuteron magnetic moment is 
underlined in table [6l since this observable is used to constrain one of the two EEC’s 
in the isoscalar current at N3LO. As already noted, loop corrections at N4LO in the 
charge operator (see figure 0]) are isovector and do not contribute to the charge radius 
and quadrupole moment. Among the N2LO and N3LO contributions to Q{d), the 
dominant one is from the two-body charge operator associated with one-pion exchange, 
panel (c) of figure 01 Lastly, the cutoff dependence represented in table |6] by the numbers 
in parentheses, is modest for these observables. 

The CST numbers for the magnetic and quadruple moments [16711174) reported in 
the table are the recent results for model WJC-2. These include the new interaction 
currents needed to conserve current that arise from the momentum dependent off-shell 
couplings that depend on the parameter Ug. As discussed in section 13.3.31 these off- 
shell couplings are needed to give a high precision fit to the np data, and also predict 
the correct three-body binding energy. The charge radius is from the older 1995 
calculation |173) and will be replaced once the new results of models WJC-1 and WJC-2 
are obtained. Note that both the magnetic and quadrupole moments are predicted by 
this calculation, and agree with experiment to about 1%. 

4-2.2. ^He and nuclei The charge and magnetic form factors obtained using the 
three theoretical approaches are compared to data in figure [221 The full result for the 
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Figure 22. (Color online) The predictions for the form factor from the conventional 
approach (thick solid line), the xEFT approach for cutoffs of 500 (thick dashed line) 
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conventional and xEFT approaches successfully reproduce the measured form factors 
up to Q ~ 3 fm“^; as a matter of fact, the conventional calculation agrees well with 
experiment up to the first diffraction region and beyond. On the other hand, even 
making allowance for the signihcant cutoff dependence, the xEFT results tend to 
underestimate the data beyond Q > 3 fm“^; in particular, they predict the zeros in 
both Fc and Fm at significantly lower values of Q than observed. As discussed above, 
the CST calculation is limited to the CIA-0 approximation, which omits interaction 
current contributions and does not fully treat the off-shell behavior of the vertex 
functions. Therefore, only the isoscalar magnetic contribution for model WJC-2 should 
be compared to data, and for this one observable (bottom left panel of figure [ 22 ]) the 
agreement is good. 

The charge and magnetic radii of ^He, the magnetic dipole moments of ^He and ^H, 
and their isoscalar and isovector combinations, are listed in table [71 The determination 
of the radii is affected by the unavailability of accurate data at low momentum 
transfer, and values for these radii are not given here. 

The conventional (yEFT) approach uses trinucleon wave functions obtained with 
the hyperspherical-harmonics method from the AV18/UIX [N3LO/N2LO(500) and 
N3LO/N2LO(600)] Hamiltonian. “Conventional” results include one-body, two-body, 
and three-body terms in the electromagnetic operators, while the yEFT results retain, 
beyond the LO terms, corrections up to N3LO in the current and N4LO in the charge 
operator. The yEFT results for the magnetic moments are underlined, since these 
observables have been used to £x the EEC’s in the current. Both the conventional and 
yEFT approaches lead to values quite close to the empirical ones. As is well known, 
lA (one-body or LO) predictions for the magnetic moments typically underestimate the 
data by ~ 15%. Finally, as previously noted, the CST results given here were obtained 
only in the complete impulse approximation (CIA-0); a more reliable prediction must 
await a calculation that includes isovector two-body currents induced by pion exchange. 

4-3. '^He nucleus 

Only results obtained in the conventional and yEFT approaches are available for 
'^He. The “conventional” (yEFT) results for the charge radius obtained with the 
AV18/UIX [N3LO/N2LO(500) and N3LO/N2LO(600)] are listed in tabled Corrections 
beyond lA (or LO) in the charge operator increase lA (or LO) predictions by about 
1%. The cutoff sensitivity, shown in parentheses, between N3LO/N2LO(500) and 
N3LO/N2LO(600) is at the % level for this observable. Conventional (AV18/UIX) and 
yEFT [N3LO/N2LO(500) and N3LO/N2LO(600)] predictions for the form factor are 
compared to data in figure [231 In the range up to Q ~ 8 fm“^ over which calculations 
have been carried out, there is satisfactory agreement between the conventional 
calculation and experiment. 
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Table 7. The ^He charge and magnetic radii, respectively rc(^He) and rm(^He) (both 
in fm), and trinucleon magnetic dipole moments g (in nuclear magnetons), obtained 
with the three different theoretical approaches, are compared to experimental values. 
The isoscalar (S') and isovector (V) combinations of the magnetic dipole moments are 
also listed. The CST results for A = i were obtained in CIA-0. The experimental 
^He charge radius is taken from Ref. m. while experimental errors on the g^s are 
negligible. The numbers in parentheses in the x^FT results give the cutoff dependence, 
while the xEFT results which are underlined correspond to those observables which 
have been used to fix the LEC’s in the current. 


Observable 

Conv. 

xEFT 

CST 

Exp. 

rc(^He) 

1.928 

1.962(4) 

1.879 

1.973(14) 

rm(^He) 

1.909 

1.920(7) 

2.035 

1.976(47) 


2.953 

2.979 

2.441 

2.979 

/i(3He) 

-2.125 

-2.128 

-1.648 

-2.128 


0.414 

0.426 

0.396 

0.426 


-2.539 

-2.553 

-2.044 

-2.553 


Table 8. The ^He charge radius (in fm) obtained in the conventional and xFFT 
approaches. The experimental value is taken from Ref. m- 


Radius Conv. yEFT 

Exp. 

rc(^He) 1.639 1.663(11) 

1.681(4) 


5. Conclusions 

Since the first measurements at Stanford, data on form factors of few-nucleon systems 
obtained from elastic electron scattering experiments have provided crucial benchmarks 
for testing our understanding of nuclear dynamics, i.e., of nuclear interactions and 
associated nuclear electromagnetic currents. The high Q measurements at SLAG 
and JLab have now pushed the experimental knowledge of these form factors (or of 
combinations of them) in the deuteron out to Q ~ 12 fm“^ and in ^He to Q ~ 9 
fm“^. Precision measurements at low Q from JLab and Mainz have also become 
recently available. The difficult (but successful) efforts to measure the deuteron tensor 
polarization, T 20 has made it possible to separate the three deuteron form factors out to 
Q ~ 9 fm“^. The data set is now sufficiently robust to permit the three deuteron and 
four three-body form factors to be extracted using the global analyses reported in this 
review. 

This precision data presents a challenge for nuclear theory. The three approaches 
discussed in this review assume that the nucleus consists of interacting protons and 
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Figure 23. The "‘He charge form factor as function of the momentum transfer Q (in 
fm“‘), calculated within the conventional (AV18/UIX) and x^FT (N3LO/N2LO(500) 
or N3LO/N2LO(600)) approaches is compared with experimental data. The thin and 
thick solid lines correspond to AV18/UIX results obtained by retaining, respectively, 
one-body only or both one- and two-body contributions in the charge operator; the 
dashed and dot-dashed lines bounding the lightly shaded area correspond, respectively, 
to the N3LO/N2LO(500) and N3LO/N2LO(600) results, and both include up to N3LO 
corrections in the charge operator. 


neutrons and that the effects of their internal substructure can be accounted for by 
effective forces and currents. The two-body forces are described, at long distance, by 
one-pion exchange and, at short distance, either by phenomenological terms of two- 
pion and shorter range in the conventional approach, or by two-pion (and three-pion) 
exchange and contact terms consistent with the symmetries of the strong interaction in 
yEFT, or by exchange of effective mesons in CST. These two-body forces are constrained 
to reproduce the large database of elastic NN differential and total cross sections and 
polarization observables up to energies close to the pion production threshold, and do so 
with a x^/datum ~ 1. Three-body forces, predominantly induced by two-pion exchange 
but supplemented by short-range or contact terms, are needed in the conventional and 
yEFT approaches in order to reproduce the three- and four-nucleon binding energies. In 
the CST calculations, irreducible three-body forces are not explicitly included. However, 
certain types of physical processes generated by the equations can be reinterpreted as 
three-body forces if one wants to establish connections to other approaches. A surprising 
result is that pure one-boson-exchange models (which generate no three-body forces) 
will give a precision fit to the NN data and correctly reproduce the empirical binding 
energy, provided these models inelude off-shell couplings at the boson exchange vertices. 
Since these off-shell couplings have no nonrelativistic analogue, a better way to compare 
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the CST with the other approaches is to transform the CST dynamics so that the off- 
shell couplings are replaced by an infinite tower of three-body forces with a unique 
structure guaranteed to give the same physics. These predicted three-body forces can 
then be compared to those required by the other approaches. 

Effective currents consist of a one-body component and interaction, or many-body 
currents. In all of the approaches, the one-body current is parametrized in terms of 
the measured nucleon electromagnetic form factors, but in order to ensure current 
conservation in the presence of off-shell nucleons, the CST requires an additional, 
unknown nucleon form factor, F^. In the conventional and xEFT approaches, the 
many-body components are induced primarily by the exchange of pions, but also the 
exchange of effective (and heavier) mesons, and by the excitation of intermediate low- 
energy resonances of the nucleon (like the A isobar) or due to transition mechanisms 
(like the pvry current). These many-body contributions are necessary for bringing 
theory into close agreement with experiment, especially in the case of the charge form 
factors of A =2-4 nuclei and the isovector combination of the trinucleon magnetic form 
factors (and magnetic moments). In the CST, the one-boson-exchange mechanisms will 
generate interaction currents only, but for the reasons discussed above, transformations 
of these currents that remove their off-shell couplings should produce predicted many- 
body currents that can be compared to those found in the other approaches. These 
studies are still in their infancy. 

One would have expected that probing few-body systems at very high Q might 
have revealed a new role for quark degrees-of-freedom at short distances. Except for 
their implicit role in determining the effective forces and currents (as well as nucleon 
electromagnetic form factors), the models discussed in this review do not include any 
such effects. However, the high Q {Q ^ 5 fm“^) data have unequivocally shown the 
need for a relativistic treatment of nuclear dynamics in this regime. The CST approach, 
which includes relativistic corrections to all orders, provides an excellent description of 
the deuteron data up to Q ~ 12 fm“^. While this conclusion was also evident from the 
2001 and 2002 reviews on the deuteron form factors [18111182] . it is strengthened by the 
new data and broader perspective of the present review. 

While the high Q CST calculations for the deuteron form factors seem to be a 
success (to be conhrmed by the next generation calculation in progress), the three- 
body calculations are still incomplete, and ^He calculations have yet to be attempted. 
The fully relativistic CST is probably too cumbersome to be extended much beyond its 
current scope, and further progress with this method would be aided by the development 
of some new approximation schemes that would permit it to be extended to more 
complicated systems. Pending such a development, the other approaches are the only 
way to study systems for H > 3 (and even to describe the H = 3 system fully). 

Perhaps the most significant new developments reviewed here are predictions for 
few-body form factors that are obtained from the new higher-order yEFT calculations 
of the two-body currents. Moreover, in the present review, we have reported the first 
yEFT calculation of the '^He form factor. The ability of these calculations to successfully 





Electromagnetic Structure of Few-Nucleon Ground States 


61 


predict the low Q form factors shows that our understanding can be traced directly to 
QCD, but the sensitivity of these calculations to the cutoffs, which becomes evident 
at high Q, shows that their predictive power is limited to low Q phenomena. This is 
really no surprise, since yEFT employs a perturbative expansion that explicitly breaks 
down at high Q. Perhaps this very promising technique can be extended to higher Q by 
employing some form of Hamiltonian dynamics [1] , but this is yet to be investigated. 
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